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ABSTRACT 

We present two luminous UV/optical flares from the nuclei of apparently inactive early- type galaxies 
at z = 0.37 and 0.33 that have the radiative properties of a flare from the tidal disruption of a star. In 
this paper we report the second candidate tidal disruption event discovery in the UV by the GALEX 
Deep Imaging Survey and present simultaneous optical light curves from the CFHTLS Deep Imaging 
Survey for both UV flares. The first few months of the UV/optical light curves are well fitted with 
the canonical f"^/'^ power-law decay predicted for emission from the fallback of debris from a tidally 
disrupted star. Chandra ACIS X-ray observations during the flares detect soft X-ray sources with 
Tbb — (2 — 5) X 10^ K or F > 3 and place limits on hard X-ray emission from an underlying AGN down 
to Lx(2 - 10 keV) < 10"^ ergs s~^ Blackbody fits to the UV/optical spectral energy distributions 
of the flares indicate peak flare luminosities of lO^"* — 10^^ ergs s~^. The temperature, luminosity, 
and light curves of both flares are in excellent agreement with emission from a tidally disrupted 
main-sequence star onto a central black hole of several times IO^Mq. The observed detection rate 
of our search over ~ 2.9 deg^ of GALEX Deep Imaging Survey data spanning from 2003 to 2007 is 
consistent with tidal disruption rates calculated from dynamical models, and we use these models to 
make predictions for the detection rates of the next generation of optical synoptic surveys. 
Subject headings: galaxies: nuclei — Ultraviolet: galaxies ~ X-rays: galaxies 



1. INTRODUCTION 

The tidal disruption of a star by a supermassive black 
hole (SMBH) has been proposed to be a unique probe 
for dorm ant bl ack holes lurking in the nuclei of norma l 
galaxies ([Frank fc Rees1ll976HL"idskii fc Ozernoi II1979D . 
When a star passes close enough to a central black hole 
to be torn apart by tidal forces, the stellar debris falls 
back onto the black ho le and produces a luminous accre- 
tion flare (|Rees IITgSSD . The detection of a tidal disrup- 
tion flare is the only direct signpost for dormant black 
holes in the nuclei of galaxies for which the dynamical 
signature of the SMBH cannot be measured, i.e. when 
the sphere of influence of the black hole is unresolved 
(i?.„/ = GMbii/<jI). 

Tidal disruption events are a potential new tool 
for measuring the SMBH mass function, independent 
of the scaling relations between central black hole 
mass and global host galaxy properties established 
for nearby galaxies with spatially resolved kinemat- 
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ics (iGebhardt et al. I 12000'; 'Ferrarese & Merritt ' '200Q!; 
iGraham et al. 1 120011 : [Marconi 
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scaling relations measured for nearby active galaxies 
with black ho l e ma, s ses from reverbera t ion mapping 
(IWandel et al. I Il999t iKaspi et al.1 120001 : IVestergaardI 
12002^. Tracing the coevolution of the mass of black holes 
and their host galaxies as a function of redshift is critical 
for constraining mechanisms for black hole growth and 
the formation of the host galaxy bulge, as well as the 
interplay between the m (Kauffmann & Haenhelt 1 120001 : 
iCattaneo et al. II2005I : ICroton et al. Il2006f r 

1.1. Theory 

The rate at which stars pass within the tidal disrup- 
tion radius of a central SMBH, Rt R^^irj^MBii/M^y/^, 
where rj is on the order of unity, depends on the flux 
of stars into loss cone orbits with Rp < Rt, where 
Rp is the distance of closest approach of the star's or- 
bit. Stellar dynamical in odels of the nuclei of g alaxies 
(jMagorrian fc Tremainel Il999t IWang &: Merritt I I2004D 
predict that this will occur once every 10^-10^ yr in a 
galaxy, depending on the nuclear stellar density profile 
and the mass of the central black hole. At these rates, 
stellar disruption events should contribute to the faint 
end of the X-ray luminosity function of active galax- 
ies, and may be an important mechanism for black 
hole growth in black hole s with A^bh < 2 x IO^Mq 
(|Milosavlievic et al~l |2006[ ). as well as a way to reg- 
ulate the tight correlation between Mbh and the ve- 
locity dispersion host galaxy bulge (jZhao et al. I l2002l : 
Merritt & Poon 2004). For higher mass black holes, 
however, there is a critical black hole mass above which 
Rt is smaller than the Schwarzschild radiu s (Rs), and a 
star is swallowed whole without disruption (jHills II1975I ). 
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M„it = 1.15 X lQ^[r^ /mf^MQ, where = i?*/i?o and 
= M^/Mq. 

After disruption, the steUar debris is distributed in 
total energy from — Ae < E < +Ae, where Ae — 
GMbhR*/Rp is the change in the black hole poten- 
tial across the radius of the star (|Lacv et al. 111983 ). If 
Ae ^ ef, GAh/Ri,, where et is the mean orbital bind- 
ing energy of the stellar debris, half of the mass of the 
star remains bound to the black hole, and half is ejected 
from the system. Numerical simulations show that the 
fraction of debris that is eventually accreted can be as 
little as O.IM^ due to the strong compression of the 
gas stream during the second passage through pericenter 
which ca n impart a frac tion of the gas with an escape 
velocity ([Aval et al. 11200 0). The most tightly bound gas 
has a total energy E = GM-Q-n/ra — — Ae, where Va is the 
apocentric radius, and thus ra — fc^^-R*(AfBH/^^*)^^'^, 
where fc = 1 if the star has no spin on disruption and 
fc = 3 if the star is spun-up to near break-up on disrup- 
tion. Numerical simulations predict that the star should 
be spun up via tidal interaction to a significant fraction of 
the b reak-up spin angular velocity (jAlexander fc Kumar I 
l2001f) . The start of the flare, occurs when the most 
tightly bound gas returns to pericenter in a time after 
the disruption, tjj, where to — to = 7r(rj^/(2GAfBH))^^^, 
and thus 



ta-tD = 0.nk-^/^{Rp/RT) 



3 3/2_-l, 



(1) 

where Mg = Mbh/IO^Mq ()Li et al. I [200^ . This time 
delay can be used as a direct measure of the central black 
hole mass, since Mbh oc k^{to — t^Y . 

The accretion rate from the fallback of the most 
tightly bound debris is predicted to result in a flare 
which peaks close to the Eddington luminosity, with a 
blackbody spectrum with T^ff ~ {LEdd/4:TTaR^y/^ = 

2.5 X lO^M^^^^r^^^ ^ K, which peaks in the extreme- 
UV (jUlmer I I1999D . The rate at which the re- 
maining debris falls-back to pericenter is dM/dt = 
{dM/de){de/dra){dra/dt) (|Phinnev 111981 . If the mass 
is equally distributed in binding energy, which is shown 
to be a r easonable approximation f r om numerical sim - 
ulations (lEvans fc Kochanek 1 119891 : lAval et al. I l2000f ). 
i.e. dM/de is constant, then dM/dt cx (GMbh)^/^^"^/^- 
During the fallback phase, the rate of gas supplied to the 
black hole follows this (t — tu)"^^'^ power-law decay, and 
determines the luminosity of the flare over the following 
months and years. 

Thermal emission is predicted to be produced during 
fallback from stream-stream collisions of the disruption 
debris. Relativistic precession of the debris orbits causes 
the ii igoing and outgoing g as streams to colhde (Rccs 
ll988HCannizzo et al. |[l990l : iKim. Par, fc Lee |[l999ll . re- 
sulting in strong shocks that quickly circularize the de- 
bris. The timescale for the debris to shock and cir- 
cularize and form an accretion torus should be on the 
order of a few times the minimum period of th e de- 
bris, given in equation 1 (jCannizzo et al. I ll99C)f) . At 
later times, the stellar debris is expected to spread out 
and form a thin disk that accretes slowly via viscous 



processes, with a mass accretion rate with a shallower 
t — tD)~^ power-law decay that approaches n = 1.2 
Cannizzo et al. |[l99Cl( ). It has also been proposed that 
an extended envelope will form around the disk from the 
expelled debris and reprocess the radiation with an effec- 
tive temperature oiT^ff » 1.3 x 10-^(M bh/10^M(^) ^"^ K 
(|Ulmer. Paczvnski. fc Goodman |[T998h . 

1.2. Previous Candidates 

Although the theoretical work described above makes 
predictions for the rate, luminosity, temperature, and de- 
cay of tidal disruption events, there is limited observa- 
tional evidence to test these predictions. The most con- 
vincing cases for a stellar disruption event occur from 
host galaxies with no evidence of an active galactic nu- 
cleus (AGN) for which an upward fluctuation in the ac- 
cretion rate could also explain a luminous UV/X-ray 
flare. A UV flare with L ~ 10"^^ ergs s^^ from the nu- 
cleus of the elliptical galaxy NGC 4552 was proposed to 
be the result of the accretion of a tidally stripped stellar 
atmosp here (iRenzini et al. Ill995f) by the mini A GN at its 
center (|Cappellari et al. lll999HXu et al. l[200l . A more 
exotic candidate for a tidal disruption event was discov- 
ered in the central region of a rich cluster: a transient 
optical source whose colors and light curve were best ex- 
plained as either a peculiar undcrluminous Type la su- 
pernova (SN la) at the redshift of the cluste r {z k, 0.25), 
or a l ensed tidal disruption event at z k. 3.3 (jStern et al. I 
[200l . 

The ROSAT All Sky Survey (jVoges et al. l[T999h con- 
ducted in 1990-1991 sampled hundreds of thousands of 
galaxies in the soft X-ray band, and detected luminous 
{lO^'^-lQ'^'^ ergs s-i), soft [Tbb = (6-12) x 10^ K and F = 
3 — 5] X-ray flares from several galaxies with no previous 
evide n ce for AGN activity (jBadc, Komossa, & Dahlcin 
1996'; 'Ko mossa fc Bade I Il9! 
r999.: Grupe. Thomas, fc Leighl 



1999t IKomossa & Greincr 
xJ [1999; Greincr et a l" 



2000t ) which were explained as tidal disruption events 



(see review by Komossa et al. 2002), an d with a flare 
rate of « 1 x 10~^ yr~^ per galaxy (jPonlev et al. I 
l2002l ). which is in agreement with the theoreti- 
cal stellar disruption rate. A decade later, follow- 
up Chandra and XMM-Newton observations of three 
of the galaxies demonstrated that they had faded 
by factors of 240 - 6000. Ahhough this dra- 
matic fading is consistent with the (t — tn)~^f^ de- 
cay of a tidal disruption flare ([Komossa et al. I 2004 



iHalpern. Gezari. fc Komossa 1 120041: iLi et al 1 120021 ) . the 

follow-up observations were obtained too long after the 
peak of the flares (5-10 yr) to measure the shape of the 
light curves. Follow-up Hubble Space Telescope (HST) 
Space Telescope Imaging Spectrograph (STIS) narrow- 
slit spectroscopy confirmed two of the galaxies as inac- 
tive, qualifying them as the most convinc ing hosts of a 
tidal disruption event (jGezari et al. l[20M l. 

The soft X-ray outburst from NGC 5905, which has the 
best sampled light curve of the ROSAT candidates, was 
fitted with a power-law decay (Komossa & Bade 

1999), and the energetics of the flare was successfully 
modeled as emission fr om the tidal dis ruption of a brown 
dwarf or giant planet ()Li et al. II2002I ). Alternative sce- 
narios for the source of the flare in NGC 5905, such as 
an SN, gamma-ray burst, lensing event, or variable ab- 
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sorbed Seyfert nucleus, were ruled out by Komossa & 
Bade (1999). NGC 5905 was subsequently found to have 
a persistent low-luminosity Seyfert nucleus in its STIS 
spectrum (Gezari et al. 2003). However, the amplitude 
of its X-ray variability (> 1000) is better attributed to 
a tidal disruption event than extreme variability of its 
low-luminosity active nucleus. 

It should be noted that the presence of a low- 
luminosity AGN and the occurrence of a tidal disruption 
event are not mutually exclusive. It is possible that the 
presence of an accretion disk may even enhance the tidal 
disruption rate due to interactions between a star and 
disk that decrease its angular momentum and bring it 
into a loss-cone orbit (ISver et al. |[l991l : lArmitage et al. I 
11996': Don lev et al. Il2002f) . Optical spectroscopic surveys 
reveal that 20 — 40% of nearby galax ies host a Seyfert nu- 
cleus (|Ho et al. Ill997t iMiller et al. | [2003). Low levels of 
AGN activity may be a ubiquitous property of galaxies 
with a central SMBH. However, the presence of an active 
nucleus does make the interpretation of a UV/X-ray flare 
subject to careful analysis. 

Two new candidates from the X MM-Newton Slew Sur - 
vey (XMMSLl) were reported by lEsquej et aP ()2006l ) 
to be optically non-active galaxies that demonstrated a 
large-amplitude (80-90) increase in soft (0.1-2 keV) X- 
ray luminosity compared to their previous ROSAT PSPC 
All-Sky Survey upper-limits, with XMMSLl luminosities 
of W^^ - 10^3 ergs s-i. 

With the exception of NGC 5905, all of the X-ray tidal 
disruption event candidates described above were simply 
"off" -"on" detections, with no detailed light-curve infor- 
mation. We initiated a program to take advantage of the 
UV sensitivity, large volume, and temporal sampling of 
the Galaxy Evolution Explorer (GALEX) Deep Imaging 
Survey (DIS) to search for stellar disruptions in the nu- 
clei of galaxies over a large range of redshifts and attempt 
to measure the detailed properties of a larger sample of 
tidal disruption events. 

In lGezari etaTI (|2006D we reported our first detection 
of a candidate tidal disruption flare by GALEX: a UV 
flare from an inactive early-type galaxy ai z — 0.3698 
with a spectral energy distribution (SED) from the opti- 
cal through the X-rays fitted by a blackbody temperature 
of (1 — 5) X 10^ K, a blackbody radius of « 10^^ cm, a 
bolometric luminosity of lO''^ — 10^^ ergs s"^, and with a 
UV light curve well-fitted by a {t — to)'^^"^ power-law, in 
excellent agreement with theoretical predictions for emis- 
sion from the tidal disruption of a star. In this paper we 
strengthen this interpretation even further by extending 
the UV light curve of the flare with another epoch of 
GALEX observations and by presenting its simultane- 
ous optical light curve from the CFHTLS Deep Imaging 
Survey. We also report a new candidate tidal disrup- 
tion flare detection from an inactive early-type galaxy 
at z = 0.326 by GALEX and present its simultaneous 
optical light curve from the CFHTLS. The broadband 
properties of the two flares, from the optical through the 
X-rays, are in excellent agreement with the predictions of 
basic tidal disruption theory and are both the most de- 
tailed observations of candidate stellar disruption events 
to date. 

The paper is organized as follows. In §2 we describe our 
techniques for selecting tidal disruption flare candidates. 



describe our sample of sources classified as quasars, and 
report the detection in the far-ultraviolet (FUV) and 
near-ultraviolet (NUV) of an SN Iln at z = 0.189 discov- 
ered by the CFHT Supernova Legacy Survey (SNLS); in 
§3 we describe our tidal disruption flare candidates in the 
UV and optical; in §4 we interpret the properties of the 
flares in the context of tidal disruption theory; in §5 we 
use the observed properties of the flares to measure the 
efficiency of our search and calculate the expected tidal 
disruption flare detection rate from theoretical models; 
and in §6 we predict the detection rate of tidal disrup- 
tion events by the next generation of optical synoptic 
surveys. Throughout this paper, calculations are made 
using Wilkinson Microwa ve Anisotropy Probe ( WMAP) 
cosmological parameters (|Bennett et al. I [20031 ): Hq — 
71 km s~i Mpc-\ = 0.27, and = 0.73. Rest- 
frame magnitudes are estim ated using the IDL r outine 
KCORRECT version 4.1.4 ()Blanton et al. |[2003[) . and 
when necessary, CFHTLS magnitudes are converted to 
the SDSS sys tem us ing the offset corrections measured 
by the SNLS (|Remau lt 2007). 

2. SELECTION OF CANDIDATES 

Our aim is to find UV flares from the nuclei of normal, 
quiescent galaxies, for which the flare is a one-time ac- 
cretion event for an otherwise dormant or undetectable 
central black hole. Such flares are best explained as the 
result of the disruption and accretion of a star. We select 
our tidal disruption flare candidates from variable UV 
sources in the GALEX Deep Imaging Survey (DIS) that 
have optical matches in the CFHTLS Deep Imaging Sur- 
vey with the optical colors and morphology of a galaxy. 
We use available X-ray catalog data to rule out galaxy 
hosts with hard X-ray emission from an AGN and obtain 
optical spectroscopy to eliminate galaxies with emission 
lines that manifest Seyfert activity in their nuclei. We 
also match the variable UV sources to the list of vari- 
able optical sources from the CFHT SNLS to check for 
anti-coincidence with SN triggers. 

The DIS covers 80 deg^ of sky in the FUV (Ae// = 
1539A) and NUV (A^// = 2316A) with a total expo- 
sure time of 30-150 ks, which is accumulated in vis- 
its during ~ 1.5 ks eclipses, when the satellite's 98.6 
minu t e orbit is in the shadow of the Earth i jMartin et al. I 
120051: iMorrissev et al. Il2005l I2007D . Due to target visi- 
bility and mission planning constraints, some 1.2 deg^ 
DIS flelds are observed over a baseline of 2 — 4 yr to 
complete the total exposure time. This large range in 
cadence of the observations allows us to probe variabil- 
ity on timescales from hours to years. Here we present 
four DIS fields that overlap with the optical Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) Deep 
Imaging Survey^: XMM/LSS (CFHTLS Dl), COSMOS 
(CFHTLS D2), GROTH (CFHTLS D3), and CFHTLS 
D4. 

In order to optimize our sensitivity to UV flares that 
decay on the timescale of months to years, we co-added 
the visits of each field into yearly epochs using the 
GALEX pipeline. Table [1] gives the coordinates and ex- 
posure times of the yearly co-added images in the FUV 
and NUV for each field from 2003 through 2007. To 
avoid the edges of the images that have degraded as- 

* Described m |http://www.cfht.hawaii.edu/Science/CFHLS| 
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trometry and photometry, as well as edge reflections, 
we only use sources with r < 0?55 from the center of 
the field. We search for variable sources in the FUV 
co-added images, instead of the NUV, (1) in order to 
avoid blending between sources, which is a larger is- 
sue in the NUV; and (2) for better contrast between 
the variable source and the host galaxy, since galaxies 
are fainter in the FUV. We use aperture magnitudes 



6") from the pipehne-g enerated SExtractor cat- 



alogs (jBertin fc Arnouts1ll996D and match them to the 
closest source within a 2" radius in the other co-added 
images, a radius that is less than half of the angular res- 
olution in the FUV (FWHM= 4.5"). We only compare 
co-added images of a field with texp > 5 ks, and measure 
the standard deviation of sources that are detected in all 
the co-added images from the weighted mean magnitude, 

o-((m)), where (m) = J2iLi ^ / J2iLi ^ is the num- 
ber of co-added images, rrii are the matched magnitudes, 
and CTi are the Poisson errors. 



2.5 1 ifrn + BskyNpix)Tc 



exp 



fm In (10) 



exp 



(2) 

where /,„ is the aperture flux in counts s^^, i?sky is the 
sky background at the position of the source in counts 
s~^ pixel"^ converted from the pipeline value for the sky 
background in counts arcsec"^ (1 pixel = 1.5^ arcsec^), 
and Npix = 16n for a 6" radius aperture (1.5" pixel"^). 
The measured standard deviation of matched magni- 
tudes between co-added images is dominated by photo- 
metric errors. For comparison, we also determine aphot 
for matched sources. 



^=l <^f 

N 



<Jphoti{m)) = 



(3) 

and the mean (ip/iot is measured as a function of magni- 
tude in bins of magnitude with at least 100 points. The 
curves are in good agreement with each other; however, 
to be conservative, we use the measured a in case there 
are systematic errors between co-added images not ac- 
counted for by aphot- 

We use a 5 cr cutoff to select intrinsically variable 
sources. The source counts start to fall off after the lim- 
iting magnitude of the shallowest coadd, so we do not 
use a measured for magnitudes fainter than mum- We 
define the limiting aperture magnitude of an exposure by 
the 5(7 detection limit for a point source: 



mii 



= -2.5 X log 5 



T, 



zp. 



exp 



(4) 

where Bsky = 3 x 10^"* cts/sec/pix is the typical sky 
background in the FUV, Bgky = 3 x 10^^ cts/sec/pix 



in the NUV, zp=18.82 in the FUV, and zp=20.08 in the 
NUV. For this step in the analysis, we do not include 
an aperture correction for the fraction of total energy 
enclosed in the aperture. 

A large, extended galaxy can be shredded by SExtrac- 
tor into multiple sources. In this case, a source from a 
shredded galaxy may not have a match in another co- 
added image if the galaxy has not been "shredded" in 
the same way in the other coadd. To avoid this problem, 
when a source in one co-added image (Ci ) has no match 
in another co-added image (C2), we place an aperture on 
C2 to measure the aperture magnitude. If the measured 
magnitude in C2 is brighter than mum of C2 then we 
assume that the source in Ci has been shredded from an 
extended galaxy, and that the lack of a source match in 
C2 is not due to variability. 

The 5cr UV variable sources are then matched with 
the TERAPIX CFHTLS-T0003 Deep Survey photomet- 
ric catalog in u [X^ff = 3743 A), g {Xeff = 4872 A), 
r {X^ff = 6282 A), i {Xeff = 7776 A), and z {Xeff 
11702 A) and the optical colors and morphology of the 
matched sources are used to identify the hosts of the 
flares. We use the closest match within a 3" radius, 
which has been determined from other GALEX studies 
to be the optimum radius to match with optical catalogs 
of similar depth, with only 1% of UV sources with mul- 
tiple optical matches within that radius ( Agii eros et al. I 
[20051) . Tables [2H5] list the FUV and NUV AB magni- 
tudes in each co-added image of the UV variable sources 
that have optical matches, as well as the factor times 
(T(m) that the source varies. The magnitude errors given 
are la statistical errors, and do not take into account 
the 2% and 5% zeropoint plus flat fielding uncertainty 
in th e NUV and FUV, respectively (Morrissey ct aQ 
|2007( ). The aperture magnitudes have been corrected for 
the fraction of energy enclosed in a 6" radius aperture, 
mtot = map + 2.51og(/ap), where fap is 0.95 and 0.86 for 
the FUV and NUV, respectively, measured from a com- 
posite point spread function (PSF) constructed in a DIS 
coadd. Sources are flagged as flares (F) in Tables [2H5] 
if they have a constant flux within the errors before the 
FUV flare, the flare fades monotonically thereafter, and 
it docs not fade fainter than its flux before the flare. If 
a monotonic FUV flare is undetected before the start of 
the flare, it is flagged in the tables as a transient flare 
(T). 

We designate optically resolved sources as those with 
a half-light radius (ri/2) > 0"6. One source (D3-15) is 
mis-classified as resolved by the optical catalog due to 
its bright magnitude (m^ < 15) which produces bright 
diffraction spikes, and this has been corrected in the ta- 
ble. Optically unresolved sources with g — r < 0.6 and 
u — g < 1 are classified as quasars, unresolved sources 
with u — g > 1.75 01 g — r < 0.6 and u ~ g > 1 are clas- 
sified as stars, and all resolved sources are classified as 
galaxies. X-ray catalog data from XMM- Newton for the 
XMM/LSS fiddJpiiaDDctti.et_alJ[200l) and the COS- 
MOS field (iManieri et al. ii2007l) and f rom Chandra for 
the GROTH field (|Nandra et al. II2005D are used to clas- 
sify galaxies with a hard X-ray point-source detection as 
AGNs. 

In addition, we match the sources to the list of spec- 
troscopically confirmed quasars and active nuclei from 
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Fig. 1. — GALEX variable sources {open circles) with CFHTLS 
optical matches {red filled circles). The blue circle shows 0?55 
deg radius of GALEX, and the green square shows 1 deg square 
diameter of CFHTLS. The area of overlap between the surveys is 
0.806 deg2 (Dl), 0.565 deg^ {D2), 0.802 deg^ (D3), and 0.709 deg^ 
(D4), for a total overlap of 2.882 deg^. The area of overlapping 
X-ray catalog data is shown shaded in magenta, and the perimeter 
of spectroscopic data from the DEEP2 Survey in the GROTH field 
is shown with a dashed black line. 



iVeron-Cettv fc VeronI (pOOl and the COSMOS Survey 
(jTrump et al. 1 120071 ). We also match our sources with 
the list of quasars identified from GALEX UV grism 
spectra in the COSMOS and GROTH fields (T. SmaU 
2007, private communication). Figure [T] shows the 
UV variable sources and their optical matches in four 
GALEX DIS fields and the overlap between GALEX and 
CFHTLS for each field, resulting in a total overlapping 
area of 2.882 deg^. Also shown is the overlap of the avail- 
able X-ray catalog data and the spectroscopic DEEP2 
survey in the GROTH field. 

The CFHTLS fields are observed up to 5 times a month 
during the seasonal visibility of each field. Real-time dif- 
ference imaging is performed on the images to produce 
a list of optically variable sources from which the CFHT 
SNLS selects SN candidates for foll ow-up spectroscopy 
(|Astier et al. 1 120061: iSullivan et aTl |2006| ). We match 
the publicly available variable optical source list (which 
includes optical variability associated with quasars and 
AGNs) with our UV variable sources. Contamination by 
SN la and gamma-ray burst afterglows is unlikely since 
they are intrinsica l ly faint in the UV (IPanagial l2003t 
iBrown et all 120051 iM^szaros fc Re"es1ll997D. How ever. 
SN H are UV brigh t at early times (jBrown et al. I [20071 : 
llmmler et aT~l |2007[ ) and are a potential contaminant, 
and so we check for anti-coincidence with the SNLS SN 
triggers. 

Figure [5] shows the optical color-color diagram of the 
UV variable source hosts and the regions of the dia- 
gram dominated by template CFHTLS colors for main- 
sequence stars, quasars, and early-type elliptical and 
spiral galaxies from the photometric redshift code Le 



I 0.5 




Fig. 2. — Optical color-color diagram of CFHTLS matches 
to GALEX UV variable sources. Resolved optical hosts are plot- 
ted with open circles, and unresolved optical hosts are plotted 
with filled circles. Sources detected by CFHT SNLS as variable 
sources in the optical are plotted with stars. Sources with an X- 
ray catalog detection are plotted with crosses, and sources that 
are spectroscopically confirmed AGNs or quasars are plotted with 
squares. Dashed lines delimit the regions of the diagram dominated 
by quasars, early-type galaxies, and main-sequence stars. 

PHARE.^ Our UV variability selection detects mostly 
quasar hosts, many of which are also detected as hard 
X-ray sources and as optically variable sources. Four of 
the optically unresolved hosts have optical colors of a 
main-sequence star; however, three of these sources (Dl- 
4, D2-6, D3-1) have a > T.'O separation from the UV 
source, which is larger than the positional uncertainty of 
GALEX and may be a false match. Our analysis does 
not detect any RR Lyrae stars, whic h are located below 
the main sequence with u — g ~1.15 (|Sesar et al. I[2007i ). 
or M-dwarf flare stars, lo cated at the top of th e main se- 
quence with u — (7 ~2.6 (jHawlev et al. Il2002f) . since our 
deep co-added images average out their variability sig- 
nal which occurs on much s horter ti mescales of days and 
minutes, respectively (Welsh et al. Il2005> ^2006) . 

Tables [BHSl give the optical AB magnitudes, half-hght 
radii, match separation of the flare hosts, and their clas- 
sification described above. Sources are flagged in the 
tables with an 'o' if they are optically variable, with an 
'x' if they are detected in an X-ray catalog, with an 's' if 
they are confirmed AGNs or quasars from optical spec- 
troscopy, and with a 'g' if they have a quasar spectrum 
in the UV. The redshift of the host galaxy is also listed 
in the table for those with spectroscopy. The UV vari- 
ability of D2-21 is due to an artifact in the 2005 COS- 
MOS coadd, and so it is labeled in Table [7] as an arti- 
fact (ART). Due to errors in the centroid of the faint 
UV source D3-14, it was erroneously not matched with 
its counterpart in the co-added images after 2003 which 
had centroids that are offset by 2.8". This source is mis- 
takenly detected as a variable source and is labeled as an 
error (ERR) in Table [3 

We define our tidal disruption event candidates, listed 
in Table [TOl as monotonic or transient UV flares from 
optically resolved galaxy hosts that are not classified as 
an AGN (no hard X-ray detection and no existing AGN 

^ Arnouts & Ilbert, [http: //www.oamp.fr/people/arnouts/LE_PHARE.html| 
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optical spectrum) . We follow up these galaxy hosts with 
optical spectroscopy, to search for signs of Seyfert activ- 
ity in their optical spectra in the form of emission lines. 
If a galaxy does not have a Seyfert nucleus, and thus no 
obvious signs of persistent accretion activity by its cen- 
tral black hole, the source of the UV flare is interpreted 
in the context of the tidal disruption of a star by a central 
black hole that is otherwise starved of gas and dormant. 

2.1. Quasars Selected from UV Variability and Optical 
Colors and Morphology 

As a biproduct of this study we have a list of 54 quasars 
classified from UV variability and optical colors and mor- 
phology alone. Only 11% of these quasars are detected 
as optically variable sources by the SNLS database. This 
discrepancy should not be due to difficulties in detecting 
variable sources in the nuclei of galaxies with image sub- 
traction, since these sources are by definition optically 
unresolved. The SNLS optical variability search is opti- 
mized for detecting SNe, which vary on the timescale of 
weeks, and within the seasonal window of observability 
of the fields. Longer timescale variability (months-years) 
may not be picked up by their selection techniques. Our 
study demonstrates that UV variability in combination 
with optical colors and morphology is an efficient method 
for detecting quasars over a large field of view, without 
the expense of X-ray surveys and optical spectroscopy. In 
Figure[3]we plot the low-state and high-state FUV— NUV 
color for all the variable sources, with quasars and AGNs 
labeled. Within the errors, all of the quasars and AGNs 
have a bluer FUV— NUV color during their high state 
in flux. This is in agr eement with the findings of opti- 
cal variability studies (|Giveon et al. I Il999l : iGeha et al. I 
[2003). This may also contribute to the better efficiency 
in detecting quasars with variability in the UV than in 
the optical. 

2.2. Detection of a Supernova Type Iln 

Candidate D4-7 is an SN candidate that was discov- 
ered by the CFHT SNLS on 2003 August 21 (SNLS ID: 
03D4ck), and was spectroscopically confirmed by Gem- 
ini on 2003 November 20 as an SN Iln at z = 0.189 
(Howell et al. 2005). GALEX detects the SN in the 
2003 co-added image of CFHTLS D4 constructed from 
visits from 2003 August 30 to 2003 September 21, with 
FUV = 23.43 ± 0.08 mag and NUV = 22.50 ± 0.04 mag. 
This is the first UV detection of an SN Iln, and one of 
onl y 5 other published UV fight curves of Type II SNe 
(see lBrown et al.1l2007l ). In a separate paper (L. Dessart, 
in preparation) we present an analysis of the detailed 
UV light curve, along with the CFHT SNLS optical light 
curve and spectra. 

2.3. Optical Spectroscopy of Candidates 

Table [TT] gives the log of spectroscopic observations 
of the host galaxies of 8 out of 10 of our tidal disrup- 
tion flare candidates, as well as their redshifts, ([O HI] 
A5007)/(H/3) narrow-line ratio, and [O HI] A5007 lu- 
minosity. We deflne a spectrum as a Seyfert 2 (Sy2) 
in the table if it has a diagn ostic narrow-line ra- 
tio (fO HI] A5007)/(H/3) > 3 (iBaldwin et al I fl98ll : 




Fig. 3.— FUV-NUV color of GALEX variable sources in their 
low state and high state of flux. Sources classified as quasars are 
plotted with an open circle, and sources classified as AGNs are 
plotted with an open square. All of the quasars and AGNs lie below 
the line of constant FUV— NUV color within the errors, indicating 
a bluer spectrum during high states of flux. 




iVeilleux fc Osterbrock II 19871: iKewlev et al. Il2001h and as 

a Seyfert 1 (Syl) if it shows broad H/3 emission. The 
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Fig. 4.— MDM 2.4 m CCDS spectrum of the galaxy host of D3-3 
that shows Seyfert-like narrow emission lines. This galaxy is also 
detected in an archival Chandra X-ray image with Lx ~ 9-3 X 10^^ 
ergs s"-*-. 

spectra are shown in Figures HHZl Dl-1, D2-9, and D3- 
3 have narrow-line ratios indicative of a Seyfert spec- 
trum . D3-3 is a Seyfert 2 galaxy at z — 0.355 that is 
also detected in an archival Chandra ACIS image in 2002 
April 10 indicating a hard X-ray luminosity of Lx{2-10 
keV) ~ 9.3 X lO'*^ ergs s"!. Three of the galaxies, D2- 
11, D3-8, and D3-21, have 2 <([0 HI] A5007)/(H/3)< 3, 
which places them in a region of the diagnostic diagram 
that needs other line ratios to determine the class of the 
galaxy. However, D2-9 and D3-21 have a broad H/3 line, 
and so we classify them as Seyfert 1 galaxies. Dl-9 and 
D3-13 have spectra with only stellar absorption lines, and 
are classified as early-type galaxies. This leaves D3-8 as 
the only galaxy without a classification from its optical 
spectrum. 

Table [TU] lists all of our candidates, the type of FUV 
fiare (monotonic or transient), whether they were iden- 
tified as optically variable in the SNLS database, their 
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Fig. 5. — Keck LRIS spectra of galaxy hosts of tidal dis- 
ruption flare candidates shown plotted at their rest wavelength. 
Hatch marks show the location of the strong O2 telluric B- and 
A-band absorption features that are not corrected for in the spec- 
tra. Dotted lines show the wavelengths of [O II]A3727, H/3, and 
[O III]AA4959,5007. All three galaxies, Dl-1, D2-9, and D2-11 have 
narrow H/3/[0 III]A5007 ratios and/or broad H/3 emission lines that 
classify them as Seyfert galaxies. 



1,0 

J 0-5 
'e 0.0 

u 

"» ^ 
S 2 

t 1 
o 

5 1,0 

0.5 

0.0 



4000 5000 6000 7000 

Rest Wovelength (A) 



Fig. 6. — Keck LRIS spectra of galaxy hosts of tidal dis- 
ruption flare candidates shown plotted at their rest wavelength. 
Hatch marks show the location of the strong O2 telluric B- and 
A-band absorption features that are not corrected for in the spec- 
tra. Dotted lines show the wavelengths of [O II]A3727, H/3, and 
[O III]AA4959,5007, [O I]A6300, Ho, and [S II]A6716. D3-8 has a 
narrow H/3/[0 HIJA5007 ratio that requires other narrow-line ra- 
tios to classify it. D3-13 shows no sign of emission lines consistent 
with its classification as an inactive early-type galaxy by its AEGIS 
Keck DEIMOS spectrum. D3-21 has a broad H/3 emission line that 
classifies it as a Seyfert 1 galaxy. 




5000 6000 7000 8000 

Wavelength (A) 

Fig. 7.— VLT FORSl spectrum of the galaxy host of Dl-9 
obtained on 2006 September 26 smoothed by 4 pixels (~10 A). The 
spectrum shows strong stellar absorption lines typical of an early- 
type galax y. The best-fitting Bruz ual-Charlot early-type galaxy 
template fBruzual fc Chariot 1120031) with an old stellar population 
and no star formation is shown in green scaled to the continuum 
of the spectrum. 

photo-z and b est-fitting gala xy tem plate from CFHTLS 
ijllbert et al l [2006; ,Lc Brun1 l2007[ ). as well as the spec- 
troscopic redshifts and classification for candidates with 
spectra. The photo-z template fitting accurately clas- 
sifies the early-type galaxies and emission-line galaxies, 
with the exception of D3-3, whose weak optical emission 
lines are not detected in the photo-z fitting. Dl-9 (Fig- 
ure [7]) and D3-13 (Figure O show no emission lines in 
their spectra, and strong stellar absorption features and 
thus are classified as inactive, early-type galaxies. We 
place upper limits on their [O IIIJA5007 luminosities of 
< 5 X 10"^^ and < 4 x 10'^^ ergs s~^, respectively. 

2.4. Optical Light Curves of Candidates 

We extract optical light curves for all of our candi- 
date fiares from the CFHTLS Deep Imaging Survey. The 
optical photometry is mea sured from difference im ag- 
ing following t he method o f lAlard fc LuptorTI (|1998f l im- 
plemented by iLeguillou I ([2003^ in the g, r, i, and z 
bands. The light curves are calibrated to the Vega sys- 
tem and must be converted to AB magnitudes using 
the offsets measured for MegaPrime/MegaCam.^° Up- 
per limits are determined for a 5a point-source detec- 
tion based on the average observing conditions using 
the MegaCam Direct Imaging Exposure Time Calculator 
(DIET).^^ Figure [H] shows the simultaneous UV/optical 
light curves extracted from the CFHTLS of the candi- 
dates with emission-line galaxy hosts. For non-transient 
flares in the UV, the low-state flux was subtracted from 
the FUV and NUV fluxes. 

Only one of the candidates, D3-3, has a flare with the 
sharp rise and monotonic decay expected from a tidal 
disruption event. The color of the fiare is very red, 
with its peak optical flux in the z band, and with a 
slope from the g to the z band well fitted with a power- 
law, Fi, oc z^", with a = —0.77 ± 0.07, which is con- 
sistent with the power-law optical continuum measured 
for AGNs (-2.0 < a < -0.43; Zheng et al. 1997; Diet- 
rich et al. 2002). Its Seyfert emission lines [L([0 HI] = 
2.5 X 10"° ergs s^^], hard X-ray luminosity {Lx « lO'^^ 
ergs s~^), and power-law optical continuum are most con- 
sistent with the UV/optical flare being a result of vari- 

http: / /www. cfht.hawaii.edu/Instruments/Imaging/Megacam/generalinforn 
http: / /www. cfht.hawaii.edu/Instruments/Imaging/Megacam/dietmegacam 
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Fig. 9. — UV/optical light curve of Dl-10, a candidate without 
an optical spectrum. The low-state flux has been subtracted from 
the FUV and NUV fluxes. Arrows show upper limits. The short 
duration of the flare and its offset from the nucleus of 0.6" are most 
likely associated with an SN II. 



Fig. 8. — UV/optical light curves of candidates with emission 
lines. For non-transient UV flares, the low-state flux has been 
subtracted from the FUV and NUV fluxes. The g, r, i, and z data 
points are measured from difference imaging with CFHTLS. For 
clarity, upper-limits (plotted with arrows) are shown only for the 
FUV, NUV, and g band. 

ability of its AGN continuum. 

Our remaining candidates without optical spec- 
troscopy, Dl-10 and Dl-14, have photo-z galaxy tem- 
plates that are consistent with an emission-line galaxy, 
although the emission lines could be from star forma- 
tion or a Seyfert nucleus. Figures [H] and [TO] show the 
UV/optical light curve of the two candidates. The short 
duration of the UV/optical flare in Dl-10 of ~ 70 days, 
as well as the offset of the variable source from the nu- 
cleus of the galaxy measured by the SNLS database of 

0.6", is compatible with an SN, and not a nuclear out- 
burst. The blue color of the flare, NUV-r = -0.6±0.1, 
and its location on the edge of a galaxy indicate that it 
is most likely an SN Type II, which are detected in the 
UV and are associated with star formation in the disks 
of galaxies. The high photo-z of Dl-14 of ~ 1.4, its er- 
ratic optical light curve, and the bell shape of its light 
curve in the NUV make this an unlikely tidal disruption 
event detection, and most likely a variable AGN. The 
only remaining candidates that satisfy our criteria for a 
monotonic UV/optical flare and an inactive galaxy host 
are thus Dl-9 and D3-13. 

3. TIDAL DISRUPTION FLARE CANDIDATES 

3.1. Dl-9: Host Galaxy Properties 

Figure [7] shows the spectrum of Dl-9, which is well 
fitted by a Bruzual-Charlot early-type galaxy template 
with an old stellar population and no star formation at 
z = 0.326. The galaxy host is not detected in the FUV 
or NUV, so we can only put lower limits on its rest-frame 
NUV — r color, which is a good diagnostic for separating 
galaxies into the red (early-type, old. A/ > 3 x lO-'^^M©) 
and blue (late-type, youn g, M < 3 x lO^^M©) sequences 
(jKauffmann et al. II20031 ). The lower limit on the UV- 
optical color, NUV — r > 4, places the galaxy on the 
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Fig. 10. — UV/optical light curve of Dl-14, a candidate without 
an optical spectrum. The low-state flux has been subtracted from 
the FUV and NUV fluxes. Arrows show upper limits. 

red side of the d ividing line between the red and blue 
sequences (jWvder et al. 1 1200711 . The optical colors of 
the host, which correspond to a rest-frame m — r=2.1 and 
Mr = —20.9, are close to the peak of the red sequence 
in th e bimodal distribution of u — r colors (jBaldrv et al. I 
[200l . 

The XMM Medium Deep Survey (XMDS) did not de- 
tect an X-ray source at the position of the host galaxy on 
2002 July 25, placing an upper limit on its unabsorbed 
flux of £ 1 X 10~^^ ergs s~^ cm~^ from 0.5 to 2 keV 
and ^ 7 X 10^^^ ergs s~^ cm~^ from 2 to 10 keV, for a 
power-law index {fi, cx v~^^~^^) typical of an AGN spec- 
trum of r = 1.7 and a Galactic colum n density towards 
the s ource of Nh = 2.61 x 10^° cm^^ (|Chiappetti et al.1 
l2005f ). For a luminosity distance of dh = 1680 Mpc, 
this corresponds to upper limits on the X-ray luminos- 
ity of Lx (0.5-2 keV) < 3 x lO^^ ergs s"! and Lx (2-10 
keV) < 2 X 10^2 ergs s^^. 

3.2. Dl-9: UV/Optical Light Curve 
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Fig. 11. — GALEX images of the tidal disruption flare candidate 
Dl-9 over 5 yr of DIS eclipses co-added into yearly epochs of NUV 
and FUV observations, labeled with the mean date of each coadd. 
No source is detected in the NUV and FUV in the 2003.86 and 
2003.97 co-added images. Circles show the 6.0" radius aperture 
used to measure the magnitudes. 
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Fig. 12. — Grey scale image of the variable optical point source 
shown on top of the i-band contours of the host galaxy of Dl-9 from 
the CFHTLS Deep Imaging Survey. The variable optical source is 
coincident with the nucleus within the positional uncertainties of 
- 0'.'25. 

Figure [TT] shows the FUV and NUV yearly coadded 
images of the tidal disruption flare candidate Dl-9. We 
include a co-added image from 2007 with tcxp = 8.2 ks 
that was obtained after the selection of candidates, and 
is not listed in Table 1. The source is undetected in the 
FUV and NUV in co-added images in 2003 composed of 
visits from 2003 October 27 to November 23 and visits 
from December 20 to 22. The source appears in 2004 
December 16, indicating an amplitude of variability of 
Am ^ 2 mag. The source then decays monotonically by 
^ 1.0 mag over 3 yr to below the detection threshold in 
2007. The UV flare from Dl-9 was also detected as an 
optically variable source in the nucleus of the galaxy by 
the CFHTLS. Figure [12] shows a gray-scale image of the 
variable optical source on top of the contours of the host 
galaxy in the i band. Figure [13] shows the optical light 
curve of the flare, with the flux from the host galaxy 
subtracted off, in the CFHTLS g,r, and i bands. The 
variable optical source is undetected in the z band, with 
z > 23.2. 

The r-band photometry catches the steep rise of the 
flare from a non-detection from 2003 August 1 to 2004 
July 24 to its peak on 2004 August 14. Two r-band 
data points were removed from 2005 August 4 and 2005 
September 4 because of systematic errors introduced by a 
poor convolution kernel that under subtracted the galaxy 
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Fig. 13. — Optical light curve of the tidal disruption flare can- 
didate Dl-9 measured from difference imaging by the CFHTLS. 
Error bars show 1 a errors, and arrows show upper limits based 
on the 5 cr detection of a point source. The optical photometry 
catches the steep rise of the flare to its peak on 2004 August 14. 
The solid line shows the {t — tr))~^^^ fit to the first 5 months of the 
g-, r-, and i-band light curves, for g — r = —0.41 and r — i = —0.33, 
and the dashed line shows the best-fitting power-law decay for the 
NUV data with to fixed to 2004.0. 
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Fig. 14. — UV/optical light curve of the tidal disruption fiare 
candidate Dl-9 measured by GALEX and CFHTLS. Error bars 
show 1 a errors, and solid arrows show upper-limits. The least- 
squares fit to the first 5 months of the optical light curve of the 
fiare with a (t — t^)"" decay with n = 5/3 is shown with a thick 
solid line. The dashed line shows the best-fitting power-law index 
to the NUV data with fixed to 2004.0, with n = 1.1 ± 0.3 and 
with NUV-r = -2.0±0.2. The dot-dash line indicates the time of 
a 30 ks ToO 0.2 — 10 keV Chandra observation that detected a soft 
X-ray source. 

near the nucleus, and results in a variable source flux 
that is systematically too bright. Figure [14] shows the 
simultaneous UV/optical light curve of the flare. Due 
to the small number of counts, the FUV and NUV 1 cr 
error bars and 95% confldence upper limits are measured 
using Bayesian statistics. 

We fit the 5-, r-, and i-band light curves simultaneously 
for 5 - r = -0.41 and r - i = -0.33 with a (i - io)"" 
power-law decay for data points from the peak of the flare 
to 2005.04. For n = 5/3, the power-law expected for a 
tidal disruption event, this results in a least-squares fit 
to the time of disruption, to = 2004.00 ± 0.02 (reduced 
X^=2.8). If we allow the power-law index n to vary, we 
get to = 2003.2±0.2 and n = 3.4±0.4 (reduced x^=2.8). 
There is a large difference in tjj for the different power- 
laws, but with the same goodness-of-fit. 
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The excellent temporal sampling of the CFHTLS ob- 
servations measures the time of the peak of the flare, 
to — 2004.622, and we can determine the rest-frame 
time delay of the most tightly bound material to re- 
turn to pericenter to high accuracy. If we assume the 
more physically motivated n = 5/3 fit to the light curve, 
then {to - to) /{I + z) ^ 0.47 ± 0.02 yr, which implies 
Mbh = (1.8 ± 0.2) fc^ X IO^Mq for a solar-type star dis- 
rupted at Rt, where k depends on the spin-up of the 
star on disruption. If we use the time of disruption 
from the n = 3.4 fit, this yields a larger time delay, 
{to — i-D)/(l + z) — 1.1 ± 0.2 yr, and a much higher 
black hole mass, Mbh = (1.0 ± 0.4)fc3 x lO^M©, which 
will quickly exceed Merit for stars with some spin on dis- 
ruption (see m.l\i . 

Although the first 5 months of the flare light curve, 
from 2004.622 to 2005.04, are well described by a (i - 
to)^^^^^ power-law decay, optical photometry in the g,r 
and i bands a year after the peak of the flare and the 
UV light curve extending 3 yr after the flare indicate 
that the emission does not continue to decay with the 
same power-law, but instead with a shallower decay. We 
fit the NUV light curve with a power-law decay with 
to fixed to 2004.0 and find a best-fitting power-law in- 
dex of n = 1.1 ± 0.3 and NUV-r = -2.0 ± 0.2. This 
power-law in the UV light curve is not consistent with the 
plateau in the optical data in late 2005 and may indicate 
a change in the UV/optical colors after 2005.04 closer to 
NUV — r ^ —0.6. The dip in the g-band light curve on 
2005.04 is not seen in the other optical bands and may 
be a measurement that is below the detection threshold 
of the observations. It may also be a systematic error, 
since the errors in the CFHTLS difference imaging fluxes 
are measured from the dispersion of flux measurements 
on a given night and do not account for systematic errors 
that can occur from constructing the difference image. If 
the dip is real, then it is a significant deviation from a 
simple power-law decay. 

3.3. Dl-9: Spectral Energy Distribution 

We triggered a 30 ks Chandra ACIS-S (0.2-10 keV) 
Target of Opportunity (ToO) X-ray observation of the 
flare on 2006 November 12 and detected four soft photons 
with energies between 0.2 and 0.4 keV (shown in Figure 
[T5|) . The background counts with energies <1 keV have a 
mean of 0.034 counts arcsec"^. For a Poisson distribution 
with a mean count rate within rap — 1.968" of 0.4 counts, 
the probability of having 4 counts from the background 
is < 0.07%, and so we consider this source a statistically 
significant detection. Our observation also puts an upper 
limit on the unabsorbed hard X-ray luminosity, L(2-10 
keV) < 6 X 10^^ ergs s~-^, for an AGN power-law F — 1.7 
and Nh = 2.61 x 10^" cm'^. 

We plot the SED of the flare over time from the optical 
through the X-rays in Figure [111 The UV/optical fluxes 
have been corrected for Galactic extinction using the ex- 
tinction curve of Cardelli et al. (1989). The simultaneous 
UV/optical fluxes on 2004.96 are well fitted by a black- 
body with a rest frame T^b = (5.5±1.0) x 10** K and an r- 
bandfluxof (4.7±0.6)x 10-3° ergs s^i cm^^ jjz-i. Using 



/^„(47r4)/(H-z) = (1.2±0.2)x 10^7 ergs 
s-i Hz-i in the r band and = {2-Kh/c^)v^ /{e^''/^'^ - 
l)47ri?^j, this corresponds to _R;,b — (7.6 ± 0.6) x 10^^ cm 
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Fig. 15. — Chandra ACIS-S image of the tidal disruption flare 
candidate Dl-9 on 2006 November 12. Photons with energies > 1 
keV arc plotted with open circles, and photons with energies < 1 
keV are plotted with filled circles. There is a statistically significant 
detection of four soft photons with energies between 0.2 and 0.4 
keV that are detected at the position of the UV fiare (labeled with 
a cross). The circle shows the aperture (rap=4 physical pixels) 
used to measure the count rate. 
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Fig. 16.— SED measured by CFHTLS, GALEX, and Chandra 
of the tidal disruption flare candidate Dl-9 over time. Fluxes have 
been corrected for Galactic extinction of Nh = 2.61 X 10^" cm^, 
and error bars show 2 a errors. Solid lines show the best fltting two- 
temperature blackbody spectrum, and the dashed line shows the 
individual blackbody components. The F > 3 power-law spectrum 
of the soft X-ray detection and the upper limit to the hard X-ray 
flux density at 2 keV are also shown. 

and Lboi = (jT'^AttRI^ = (3.8 ± 0.7) x lO*^ gj.gs g-i jf 
we assume that the temperature of the emission does not 
change, and we scale this blackbody to the optical fluxes 
at the flare's peak on 2004.62, this corresponds to a bolo- 
metric luminosity of (9.1 ± 0.2) x 10*"^ ergs . When 
we fit the optical fluxes at the peak with a power-law, 
Fi, cx J/", we get a = -1-1.5 ± 0.01, which is unlike the 
a < —.53 typical of AGNs and is better ascribed to the 
Rayleigh- Jeans tail of a blackbody spectrum. 

During the soft X-ray detection on 2006.88 we require 
a second higher temperature blackbody to fit the X-ray 
flux density at 0.3 keV, with Tbt = (2.7 ± 0.2) x 10^ K. 
For this blackbody spectrum, the soft X-ray count rate 
on 2006.88 corresponds to an unabsorbed (0.2-1.0 keV) 
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luminosity of (3.5 ±0.2) x IQ-^^ ergs s"^ This blackbody 
has Rbb = (6±1) X 10^2 cm and Lboi = (1.4 ±0.6) x Iff** 
ergs s~^. Due to the small number of counts, we cannot 
fit the spectrum of the X-ray emission, and so we cannot 
be sure that it is truly a soft blackbody. However, we 
can place limits on a power-law fit, Fe oc , to the 
spectrum based on the lack of a source detection above 1 
keV, to have T > 3 and L(0.2 - 1.0 keV) < 5 x 10^^ ergs 
s~^. Unfortunately, we do not have X-ray observations 
during the peak of the flare to determine if the high- 
temperature component is present and with the same ra- 
tio to the lower temperature blackbody as at later times. 

3.4. D3-13: Host Galaxy Properties 

D3-13 was r eported in iGezari et al. I (|2006) for hav- 
ing an AEGIS ijDavis et al. Il2006f) Keck DEIMOS spec- 
trum that classified the host galaxy as an inactive early- 
type galaxy at z = 0.3698. We also obtained a Keck 
LRIS spectrum of D3-13 (shown in Figure O to cover 
the wavelength region near Ha, and no Ha emission is 
detected, consistent with its previous classification as 
an inactive galaxy. The galaxy host is not detected in 
the FUV or NUV, so we can only put lower limits on 
its rest-frame NUV— r ^ 4, which places the galaxy on 
the red side of the dividing line be tween the red and 
blue sequences (|Wvder et al. II2007I ). The optical col- 
ors of the host galaxy, which correspond to a rest-frame 
u— r=2.0 and = —21.7, place it on the blue end of the 
u—r color distributi on corresponding to the red sequence 
(jBaldrv et al. Il200^ . The location of the galaxy on the 
edge of the red sequence is in agreement with the mor- 
phology of the galaxy, since the bulge-disk composition 
of its AEGIS HST ACS /-band image doe s detect a disk 
comp onent, with a bulge fraction of 0.72 (|Gezari et al. I 
[20061) . No hard X -ray source was detected in deep (190 
ks) archival AEGIS Chandra ACIS-I observations of the 
source taken over 2005 April ~ September, resulting in 
an upper limit to the unabsorbed (2 — 10 keV) X-ray lu- 
minosity of < 1.5 X lO''^ ergs s~^, for an AGN power-law 
r = 1.4 and Nr = 1.3 x 10^" cm^^ and (11 = 1970 Mpc. 

3.5. D3-13: UV/Optical Light Curve 

The UV l i ght c urve of D3-13 was presented in 
iGezari et al. I (|2006f ) , with no optical variability reported 
in CFHTLS data spanning from 2005 January to June 
with miim ~ 25. Here we extend the UV light curve with 
one more epoch of CALEX DIS data from 2007. We also 
use a deeper co-added image in 2004 with texp — 8.3 ks, 
which includes visits that have positional offsets from the 
center and thus were not included in the main coadd. We 
also present new data from CFHTLS that detect an op- 
tically variable source in the nucleus of the host galaxy 
in 2004 (Figure [T7|) . Figure [H] shows the optical light 
curve in the g, r, and i bands with the flux from the host 
galaxy subtracted off. The variable source is undetected 
in the z band, with z > 23.2. 

The r- and i-band light curves catch the flare close to 
its peak on 2004.04 and its subsequent steep decay. The 
light curves are best described with a broken power-law. 
A solid line shows the best- fitted n = 5/3 power-law 
decay to the first few months of the g-, r-, and i-band 
data from 2004.04 to 2004.35, simultaneously fitted for 
g—r = —0.11, r—i = —0.13, yielding a tight constraint on 
to = 2003.71 ± 0.01 (reduced = 3-2). When we allow 




1 0-1-2 
Arc Seconds 
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Fig. 17. — Gray-scale image of the variable optical point source 
shown on top of the i-band contours of the host galaxy of D3- 
13 from the CFHTLS Deep Imaging Survey. The variable optical 
source is coincident with the nucleus within the positional uncer- 
tainties of 0''25. 
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Fig. 18. — Optical light curve of the tidal disruption flare can- 
didate D3-13 measured from difference imaging by the CFHTLS. 
Error bars show 1 a errors, and arrows show upper limits based on 
the 5cr detection of a point source. The r- and i-band photometries 
catch the light curve close to the peak of the flare. The solid line 
shows the {t — to)~^^^ fit to the g-, r-, and i-band light curves 
for g — r = —0.11 and r — i = —0.13, and the dashed line shows 
the best-fitted power-law index for the NUV data with tj5 fixed 
to 2003.7. The dot-dash line indicates the time of the soft X-ray 
detection by Chandra. 



n to vary, the resulting fit yields tjj = 2003.686 ± 0.003 
and n = 1.77 ± 0.01 (reduced = 2.9). We include 
the sensitivity of the time of disruption to the power- law 
index in the error bars for to and use i^j = 2003.70±0.02. 

Figure fT9l shows the simultaneous FUV, NUV, and op- 
tical light curve. Due to the small number of counts, the 
FUV and NUV 1 a error bars and 95% confidence upper 
limits are measured using Bayesian statistics. The NUV 
photometry was measured with rap = 3.75" and fap — 
0.70 to avoid contamination by a pair of foreground- 
interacting galaxies detected in the NUV just south of 
the source. For the later observations in 2005 the FUV 
detector was temporarily not operational. 

The UV data trace the decay of t he flare at later times , 
which has a shallower power-law. In lGezari et al. I ()2006f ) 
we fitted the NUV light curve with an n = 5/3 power- 
law decay, and found a best fit t^) = 2003.3 ± 0.2. When 
the power-law index was allowed to vary, we found tjj — 
2002 ± 2 and n = 3 ± 2. Here we fit the optical and 
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Year 

Fig. 19. — UV/optical light curve of the tidal disruption flare 
candidate D3-13 measured by GALEX and CFHTLS. Error bars 
show 1 a errors, and solid arrows show upper-limits. The least- 
squares fit to the first 4 months of the optical light curve of the 
fiare with a(t — t£))~" decay with n = 5/3 is shown with a thick 
solid line. The best-fit power-law index for the NUV data after 
2004.04 with tn fixed to 2003.7 and NUV-r = -1.4 is shown with 
a dashed line with n = 0.82±0.03 . Tick marks indicate the time of 
archival AEGIS 100 ks 0.2 - 10 keV Chandra observations. A dot- 
dash line indicates the time of a soft X-ray detection by Chandra. 

NUV data points after 2004.35 with NUV-r = -1.4, 
the optical colors from the n — 5/3 fit, and with to 
fixed to 2003.7, and we get a best-fit power-law index of 
n = 0.82 ±0.03. 

Although the optical light curve of D3-13 puts strong 
constraints onto, we do not see the rise of the flare to its 
peak, and so we use the earliest detection of the flare in 
the r and i band on 2004 January 13 to place an upper 
limit on (to — to) /(I + z) < 0.25 yr, which corresponds 
to Mbh <5.2 xlO°fc^M0 for a solar-type star that is 
disrupted at Rt- 



I'- 



Fig. 20.— SED measured by CFHTLS, GALEX, and Chandra of 
the tidal disruption flare candidate D3-13 over time. Fluxes have 
been corrected for Galactic extinction of A'^^f = 1.30 X 10^" cm~^, 
and error bars show 2 a errors. Solid lines show the best-fltting 
two-temperature blackbody spectrum, and the dashed line shows 
the individual blackbody components. The steep slope of the soft 
X-ray detection (P = 7 zb 2) and the upper limit to the hard X-ray 
flux density at 2 keV are also shown. 

higher temperature (T2) blackbody spectrum is well con- 
strained by the soft X-ray flux detection on 2005 April 
7, with Tbb = (4.9 ± 0.2) x 10^ K, and the steep slope 
of the soft X-ray detection could be ascribed to emission 
from the Wien's tail of a blackbody. This blackbody has 
Rbb = (1.3±0.2) X 1013 cm and Ltoi = (7±l)x 10^^ ergs 
s~i. The lowest temperature for T2 that is consistent 
with the UV/optical fluxes on 2004.35 is Tbb = 1.2 x 10^ 
K and corresponds to Rbb ~ 4.7 x 10^'^ cm and Lboi > 
3.3 X 10^'' ergs a'^. 

The bolometric luminosity of the flare is dominated 
by the higher temperature blackbody component. If we 
scale the blackbody model to the optical fluxes near the 
peak of the flare on 2004.04, this indicates a peak bolo- 
metric luminosity (for T2 > 1.2 x 10^ K) of > 2.5 x lO"^ 
ergs s~i. If we fit the optical fluxes with a power-law, we 
get a — +0.5 ± 0.2, which is unlike the a < — .53 typical 
of AGNs. 



3.6. D3-13: Spectral Energy Distribution 

We plot the broadband SED of D3-13 over time from 
the optical through the X-rays in Figure[2Dl Two archival 
45 ks AEGIS Chandra ACIS-I observations were taken 
on 2005 April 6 and 7, and an extremely soft source with 
10 photons with energies between 0.3 and 0.8 keV ap- 
peared in the first ^ 30 ks of the second day of obser- 
vations, with a steep spectral slope of F = 7 ± 2 and a 
corresponding unabsorbed luminosity of L(.3 — 1 keV) 
~ 1 X lO''^ ergs s^i. During the 45 ks observation on 
2005 April 6, and during a later 100 ks observation on 
2005 September 20-23, no source was detected from 0.3 
to 2.0 keV. For a source with the same spectrum, this 
places upper limits on L(0.3 — 1 keV) of < 2 x 10"*^ and 
< 1 X lO''^ ergs s~i, indicating variability of a factor of 
10 on the timescale of 1 day. The variable soft X-ray 
source is unlikely a contaminant in the image, since the 
duration of the fiare is much longer than a cosmic-ray 
afterglow, which typically lasts on the order of hundreds 
of seconds (E. Laird 2007, private communication). 

The simultaneous UV/optical SED of the fiare on 
2004.35 requires a two temperature blackbody fit. The 
low-temperature (Ti) blackbody component has a rest- 
frame Tf,f, = (1.0 ±0.1) X 10-1 K, with Rbb = (4.7 ±0.3) x 
IQ^^ cm and Lboi = (1.6 ± 0.7) x lO'^^ grgs s^^. The 



4. INTERPRETATION 

The UV-optical colors of the candidate fiares 
(NUV-r < -1) are incompatible with SNd la and GRB 
afterglows, which are intrinsically faint in the UV, and 
their locations at the centers of the galaxies are not char- 
acteristic of SNe II which are associated with massive star 
formation in the disks of galaxies. The peak bolometric 
luminosity of the nuclear fiares (> 10**^ — 10"*^ ergs s"^) 
makes them most likely associated with an outburst from 
a central SMBH with Mbh > IO^^Mq. Chandra obser- 
vations of the host galaxies put limits on the presence 
of hard X-ray emission from a power-law AGN in the 
nuclei of the galaxies down to Lx ^ 10"*^ ergs s"^, and 
ground-based optical spectra put limits on optical emis- 
sion fines powered by an AGN to L([0 III]A5007) < 10^" 
ergs s~i. These upper limits do allow for the presence of 
a low-luminosity AGN (LLAGN) in these galaxies at the 
faint end of the lum inosity funct ion measured for nearby 
LLAGNs (Hcckm an et al. 1120051 ). In addition, the ampli- 
tude of variability in the UV in these galaxies (Am > 2) 
is comparable to that of some of the AGNs and quasars 
in our sample. However, the detailed optical light curves 
of these candidates, with their steep rise to the peak and 
subsequent power-law decay, are distinctive from AGNs, 
as are the positive optical power-law slopes of their SEDs 
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Fig. 21. — Solutions for the central black hole mass of Dl-9 as a 
function of mass for main sequence stars with a spin-up parameter 
k. The thick solid line shows the maximum possible black hole 
mass for Rt > Rsch, and the dashed line shows the maximum 
black hole mass for Rti, > Rg . 

{a > +0.5). The combination of the power-law light 
curves, blackbody SEDs, soft X-ray properties, lack of 
detectable Seyfert activity, and large luminosities make 
a tidal disruption event the most likely explanation for 
the source of the UV/optical flares. 

The change in the exponent of the power-law decay in 
the light curves of the candidates occurs (t — t£))/{l+z) 
= 0.78 and 0.47 yr after the disruption for Dl-9 and D3- 
13, respectively. These timescales are consistent with 
the timescale for circularization of the debris, which is 
expected to be on the order of the timescale of the mini- 
mum period of the orbital debris. This minimum period 
is measured for Dl-9 and D3-13 to be (to — tD)/{l+z) = 
0.47 and 0.25 yr, respectively, which would correspond to 
circularizati on in ^ 2 orbital periods. Numerical simu- 
lations from I Aval et al. I (|2000l ) also show that the accre- 
tion rate can also deviate from a (t — ^d)"^^'^ power-law 
during the fallback phase due to the expulsion of de- 
bris from compression during the second passage through 
pericenter. The broken power-law model does not de- 
scribe the plateau in the optical light curve (and dip in 
the g band) of Dl-9, or the nearly constant UV flux of 
D3-13 at late times. The emission mechanisms during a 
tidal disruption event (fallback, stream-stream collisions, 
compression, accretion) are complicated and not well un- 
derstood. We are seeing the messy details of the light 
curves of these events for the first time, however, and 
it is impressive that a simple broken power-law model 
reproduces the observations as well as it does. 

4.1. Central Black Hole Mass 

If we assume that the debris has shocked and circular- 
ized into an accretion torus or disk after the change in the 
power-law of the light curve, then the soft X-ray emission 
detected during this time may originate from the inner 
regions of the disk or torus. This interpretation would 
also account for the rapid variability observed for the 
soft X-ray source in D3-13. The radius of emission dur- 
ing the soft X-ray detections can place constraints on the 
central black hole mass, since this emission should have a 
pericenter radius that is greater than the minimum sta- 
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Fig. 22. — Solutions for the central black hole mass of D3-13 as a 
function of mass for main sequence stars with a spin-up parameter 
k. The thick solid line shows the maximum possible black hole 
mass for Rt > Rsch, and the dashed line shows the maximum 
black hole mass for ijf,6 > Rg- The horizontal arrow shows the 
minimum mass of the disrupted star derived from the minimum 
mass accreted. 

ble particle orbit for the black hole (Rms), which ranges 
from Rms = 6i?g for a black hole with no spin down to 
Rms = Rg for a maximally spinning black hole, where 
Rg — GMbb/c^- If we require that Rbb > Rms, then the 
blackbody radii of the soft X-ray components for Dl-9 
and D3-13 of 6 x 10^^ and 1.3 x lO^'^ cm, respectively, 
require that Mbh < 4 x lO^Af© and 8.7 x IO^'Mq, which 
are also consistent with Mbh < Merit for a solar-type 
star. 

There is an uncertainty in the peak bolometric lumi- 
nosity of the flares in Dl-9 and D3-13 due to the lack 
of X-ray observations close to the peak of the flare. 
However, we can estimate the minimum mass accreted 
during the fallback phase of the debris by taking the 
lower limit to the peak luminosity, and integrating un- 
der the portion of the light curve well described by a 
f-5/3 power-law to get the total energy output. The lu- 
minosity of the lower temperature component of Dl-9 
gives a lower-limit to its peak bolometric luminosity of 
9.1 X lO'*^ ergs s~^. We integrate under the portion of 
the light curve well described by a t~^^^ power-law to get 
Jl"°'-'>''Lo[{t-tD)/{to-tD)]-'/'dt = 8.0X IQSO ergss-i 

which corresponds to Etotl{e.(?) = O.OO4M0 accreted, for 
an efficiency of converting mass into energy of e = 0.1. 
This places a lower limit to the mass of the disrupted 
star to be 0.008Mq if at least half of the stellar debris 
is unbound from the system at the time of disruption. If 
we do the same exercise for D3-13, using the lower-limit 
to the peak luminosity from the lower bound to T2 of 
2.5 X 10*^ ergs s~^, and integrate from to 2004.35, this 
results in a minimum mass accreted of 0.11 Mq, which 
requires a star of at least 0.22Mq. 

In Figures [2T] and [22] we plot the solutions for Mbh 
from the rest-frame time delay between the time of the 
disruption and the peak of the flare, (to — t£))/(l + -2), as a 
function of main-sequence star mass and fc, in comparison 
to Merit 1 which is the upper limit to the black hole mass 
for which a star will be disrupted, i.e. Rt > Rsch- For 
1 < fc < 3, i.e. the star is spun-up to some fraction 
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and with an effective temperature cfiaracterized by Ed- 
dington luminosity radiation at the tidal disruption ra- 
dius, 



Teff = 2.5 X 105(Mbh/10'5Mq)1/i' K. 



(6) 



Simila rly to the calculation of iMagorrian fc Tremaine I 
(fl999ll. we derive the black hole mass function from the 
iFerguson fc Sandage I ()1991[ ) E + SO luminosity function, 



Fig. 23. — Predicted peak FUV and g magnitudes and of flares 
from 10^ — 5 X IO^Mq black holes in comparison to the detection 
limit for a 20 ks GALEX DIS exposure and a 720 s Pan-STARRS 
MDS exposure. 

of break-up, we get a range of allowable solutions for 
Mbh for Dl-9 from 8.2 x lO^M© to the upper hmit from 
Rbb > i?3 of 4 X lO^M©. The maximum allowable value 
for k in this range of black hole masses is k = 1.7. For 
D3-13 the full range of fc-values is allowed, which yields 
a large range of black hole masses from a minimum of 
2.4 X IO^A/q to the upper limit from Rtb > Rg of 8.7 x 
IO^Mq. 

We can compare these values for Mbu with the mass 
expected from the properties of the host galaxies. For 
Dl-9 the rest-frame absolute i?-band Vega magnitude is 
Mb — g — 51og((ii,) + 5 ~ K = —19.5, where is in 
Mpc and K = 0.68. This can be related to black hole 
mass using the correlat i on be tween Mbh and is (bulge) 
from IMagorrian et al.1 (|1998f ) which has a large scatter 
(~0.5 dex) that yields Mbh ~ iti^r ^ IO^Mq, which is 
compatible with solutions in Figure [211 Because of the 
high spectral resolution of the AEGIS Keck DEIMOS 
spectrum of D3-13, it is possible to measure the velocity 
dispersion (cr) of the bulge directly and use the observed 
tight correlation between Mbh and a (jOebhardt et al. I 
I2OOOD . The Mbh estimated from the velocity dispersion, 
measured to be = 120 ± 10 km s"\ is 2^1 x IO'^Mq 
(jGezari et al. I [20061 ) and is compatible for solutions in 
Figure [m with k < 2. 

5. TIDAL DISRUPTION RATE 

In this section we describe how we measure the effi- 
ciency of our search over 2.882 deg^ of DIS and simulate 
the tidal disruption flare detection rate expected from 
theor etical models . We use the tidal disruption rates 
from [Wang fc Merritt I ()2004[ ) calculated for 51 ellipti- 
cal galaxies with measured surface brightness profiles. 
We fit the galaxies with Mbh ~ IO^Mq with a power- 
law to get a black-hole-mass-dependent tidal disrup- 
tion rate A^(Mbh) ~ 1-6 x 10-'*2/r-i(MBH/10^MQ)--3". 
This gives lower rates than the analytically determined 
rate for a singular isothermal sphere with a stellar 
density profile of p oc = iV(MBH) ^ 6.5 x 

10-4y^-i(MBH/10^Mo)--25. We assume that the flares 
radiate at the Eddington luminosity, 



iEdd = 1.3 X 1044(MBH/10^AfQ)ergs s" 



N{L)dlogL 



No 



(V2^A) 



exp 



1 flogL- log Lq 
2 



dlogL, 



(7) 



where A = 0.6, Nq = 3.0 x 10"^ Mpc^^^ ^^^^ ^ 
2.6 X IO^Lq, with a factor of 2 to account for bulges in 
early- type spirals, scaled to black ho le mass by combin- 
ing th e Mbh - Mb^ige relation from , Merritt fc Ferrarese I 
([20011) and the mean g alaxy mass-to-light ratio from 
[Magorrian et al. I ()1998h which yields, Mbh = 1-6 x 



IO- 



CS) 



( j Mq. We only consider the range of 

black hole masses for which a tidal disruption flare can 
radiate close to IjEddi the 10^ — 5 x 10 ""M© mass range 
(jUlmer [[1991 ). Figure [23[ shows the predicted FUV mag- 
nitudes of flares from 10^ — 5 x lO^M© black holes in 
comparison to the detection limit for a 20 ks DIS ex- 
posure. We include an estimate of the attenuation in 
the FUV band by neut ral hydrogen ab sorption systems 
in the line of sight from [Madau I ( 1995[ ). which produces 
a non-negligible effect for sources at z > 0.6. The K- 
correction, K{z) « (-7.5)log(l-|-z), is large, and thus Ed- 
dington luminosity flares from a 5 x lO^Af© black hole 
can be detected in the FUV band out to z ^ 1.5. 

The detectability of a flare by GALEX also depends 
on the FUV magnitude of its host galaxy. We estimate 
the FUV magnitude of the host galaxy of a given black 
hole mass by converting the i3-band luminosity of the 
galaxy to a (/-band a pparent magnitude usi ng average 
g — B corrections from lFukugita et al. I (|1995[ ) measured 
for elliptical/so galaxies of —0.375 and for Sab/Sbc spi- 
ral galaxies of —0.340. We then convert the g magnitude 
to an FUV magnitude with g— FUV> 3 for elliptical/SO 
galaxies and 5— FUV> 2 for Sab/Sbc spiral galaxies. The 
number of events detected within a redshift z over a base- 
line of observations of At years is thus 



iV(< z) = AAt X 

5xlO^A/0 

N{MBli)NBH{MBli)V{z)S{MBH,z)f{MBli)dMBH, 

where A is the area in degrees, S{AIbb, z) is the selection 
function, and / is the fraction of flares that radiat e at the 
Eddin gton luminosity for a solar-type star from [uTmer I 
(|1999D . Because the peak accretion rate of the debris 

— 1/2 

is proportional to ^/bh ' central black hole mass 

increases, the fraction of flares that have a mass accretion 
rate close to the Eddington limit and will have Lp^ak ~ 
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^Edd decreases. We use / = 1 for Mbh < 4 x 10^ Mq, 
and for Mbh > 4 x lO^M© we use the log-linear function, 



/(Mbh) = -0.909 log(MBH/M0) + 7.00. 

(8) 

We estimate the selection function with a Monte Carlo 
simulation. We simulate the number of flares pre- 
dicted for each black hole bin for the redshift volume 
V{z), selected from a Poisson distribution with a mean 

(7V(Mbh,^)) = AAtN{MBu)NBH{MBH)V{z)fiMBu), 
and randomly select z and to for each flare. Val- 
ues of tQ are uniformly distributed between 2003 and 
2003-|-Ai, and the distribution of z-values is propor- 
tional to the CO- moving volume at each redshift, V{z) = 

We assume that the flares decay as 



41,253 



m{t) = Too — 2.5 log 



{t - to) 



-1 -5/3N 



- to) 



A 



FUV 



(9) 

where (to - to) = 0.11/c-3/2(AfBH/10^MQ)i/2; fc = 1,2, 
or 3; and Apuv = 8.24 x _E(i3 — T^) is the Galactic extinc- 
tion in the FUV (jWyder et al. l[2007l ) measured from the 
average Milky Way extinc tion for t he field listed in Table 
1 of Schle gel et aU ()1998f) . and the lCardelh et al. I (|1989f ) 
extinction curve. We record the magnitude of the flaring 
sources in each epoch of observations for each field with a 
co-added image with t^xp > 5 ks, and if the magnitude in 
an epoch is fainter than the limiting magnitude of the co- 
added image of the field, it is marked as undetected. The 
standard deviation is measured and compared to the 5 a 
cutoff for each field to see if it would have been selected 
by our algorithm (described in §2). We then measure 
the fraction of flares detected to determine S'(Mbh, z) for 
each field. The error in the detection rate then includes 
the Poisson error in the number of flares, plus the stan- 
dard deviation of the fraction events detected between 
10 runs of the simulation. 

In Figure [24] we plot the cumulative number of events 
that would be detected in our search within a redshift 
volume for different values of the star's spin-up param- 
eter, fc, compared to our observed cumulative detection 
of two events within z < 0.4. Each redshift volume is 
from a separate simulation, so it is encouraging that the 
cumulative distribution converges and is smoothly vary- 
ing. The k = 1 curve, corresponding to a star with no 
spin, is strongly ruled out because the slow decay of the 
flares, {to ~t£, oc k^^/^), would have resulted in an order 
of magnitude more flare detections. The fc=3 curve is in 
good agreement with our detection rate for z < 0.4, and 
the fc = 2 curve is within 1 a, but they both over pre- 
dict the number of events beyond z ^ 0.4. The detection 
rate for all values of fc levels off after z ^ 0.7, with a total 
number of detections of 17 ± 4, 6 ± 2, and 3 ± 1 for fc=l, 
2, and 3, respectively. 

The theoretical curves predict the detection of at least 
one more flare with 0.4 ^ z ^ 0.7. Figure [55] shows the 
number of events detected as a function of black hole 





2.0 



Fig. 24. — Cumulative number of events expected to be detected 
in our study within a volume out to redshift z, for different values 
of the star's spin up parameter k. The observed cumulative number 
of events detected is shown with a thick solid line. 
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Fig. 25. — Solid lines: Number of events expected to be detected 
within a volume z < 0.4, and 0.4 < z < 0.7, as a function of black 
hole mass, for A; = 3 and the fraction of flares that radiate close 
to LEdd for 1 Mq stars from Ulmer (1999). The peak sensitivity 
of our search shifts from Mbh = f-6 x IO^A/q for z < 0.4 to 
M-BH = 2.2 X IO'^Mq for OA < z < 0.7, with no sensitivity for 
Mbh < 4 X 10® Mq black holes beyond z=0.4. Dashed lines: Same 
for the fraction of flares that radiate close to Lfidd calculated for 
a more realistic Salpeter mass function from Ulmer (1999), which 
causes a decrease in number of events expected to be detected in 
the volume OA < z < 0.7. 

mass (fitted with a polynomial function for clarity) for 
the volume z < 0.4 in comparison to 0.4 < z < 0.7. 
The peak sensitivity of our search shifts from Mbh = 
1.6 X IO'^Mq for z < 0.4 to Mbh = 2.2 x lO^M© for 
0.4 < z < 0.7, with no sensitivity to black holes with 
Mbh < 4 X IO^Mq beyond z = 0.4. Thus, the deficit 
of detections in this redshift range may be the result of 
overestimating the fraction of flares from massive black 
holes that radiate at the Eddington limit. It is not sur- 
prising that the two tidal disruption flares we detected in 
our search have properties indicating black hole masses 
of a few xlO^M0, since that is where our detection sen- 
sitivity is at its peak. 
We repeated our simulation with a fraction of tidal 
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Fig. 26. — Cumulative number of events expected to be detected 
in our study within a volume out to redshift z for different values 
of the star's spin up parameter k, and the fraction of flares with an 
Eddington luminosity flare for stars with a Salpeter mass function 
between 0.3 and IMq from Ulmer (1999). The observed cumulative 
number of events detected is shown with a thick solid line. 

disruption flares that radiate at the Eddington himinos- 
ity calculated for a more realisti c Salpeter ma ss func- 
tion between 0.3 and 1.0 Mq by lUlmer I (|1999[ ). which 
gives a lower fraction for high-mass black holes than 
the curve calculated for IMq stars only (equation 8). 
We now use / = 1 for Mbh < 2.5 x lO^M©, and for 
Mbh > 2.5 X 10^ Mq we use the log-linear function, 



/(Mbh) = -0.850 log(MBH/MQ) + 6.44. 

(10) 

Figure [26] shows the resulting cumulative detection rate 
within a redshift volume, with a total number of detec- 
tions expected of 10 ± 3, 4 ± 2, and 2 ± 1 for fc=l, 2, 
and 3, respectively. The k = 2 and fc = 3 curves are now 
both consistent with the observations within 1 cr, while 
the k = I curve is still strongly ruled out. 

It is surprising that the models are in such good agree- 
ment with our detection rate, since we have not taken 
into account the effects of extinction internal to the host 
galaxy, or the inclination angle of the galaxy. Extinction 
intrinsic to the host galaxy should suppress the detection 
rate equally over the range of black hole masses, and thus 
could allow for lower values oi k, k ^ 2, to be consistent 
within Icr of our observed detection rate. This is impor- 
tant, since the analysis of the individual detections favor 
smaller values of fc. 

It appears that the values for 7V(AfBH) calculated from 
the galaxy dynamical models are in good agreement with 
our detection rates, and can be used along with the ob- 
served properties of the flares, to make detailed predic- 
tions for the detection capabilities of other surveys. 

6. DISCUSSION 

For the first time we have measured the broadband 
SED and detailed light curves of two tidal disruption flare 
candidates. The blackbody temperatures, luminosities, 
and power-law decays of the flares are in excellent agree- 
ment with that expected from the simplest "fallback" 
models for the accretion of a tidally disrupted star. 



The statistics of the detection rate of our study re- 
quire for a star to be spun up on disruption (fc ^ 2) and 
for the fraction of flares that radiate at the Eddington 
luminosity to decline for Mbh > 2.5 x 10^ A/©. How- 
ever, there is some degeneracy between the effects of in- 
ternal extinction and the spin-up parameter of the star, 
and so smaller values of k are compatible with our ob- 
servations, if we have underestimated the effects of ex- 
tinction and inclination angle of the host galaxy. The 
black hole mass-dependent tidal disruption rates calcu- 
lated from galaxy dynamical models appear to reproduce 
the observed rates, given our basic assumptions for the 
luminosity and temperature of the events. 

The detection of tidal disruption flares in the optical 
for the first time has important consequences for the de- 
tection capabilities of future large optical synoptic sur- 
veys such as Pan-STARRS and LSST. Since the the- 
oretical models for fc = 3 and /(Mbh) for a Salpeter 
mass function reproduce the observations so well, we can 
apply our analysis to predict the detection rate of the 
upcoming Pan-STARRS Medium Deep Survey (MDS), 
which will monitor an area of « 50 deg^ daily with one 
of five filters {g,r,i, z,a,nd Y), and with the same filter 
every 4 days, with a limiting magnitude in the g band 
of nig < 24.76. We measure the selection function by 
simulating flares and requiring that the flare be above 
the detection threshold of the individual exposures for 
at least 8 days, so that the power-law decays of a tidal 
disruption flare can be easily distinguished from optical 
variability from AGNs. We do not correct for Galactic 
extinction, since the MDS fields are at high Galactic lat- 
itudes. 

Figure [23] shows the predicted peak magnitude of flares 
in the g band as a function of redshift and black hole 
mass, and Figure [57] shows the detection rate by MDS as 
a function of black hole mass for fc = 3 and /(Mbh) for a 
Salpeter mass function, in a volume of z < 0.5, resulting 
in a total detection of ~ 15 events yr~^. Figure [771 also 
shows the detection rate of the Pan-STARRS Sir survey, 
which will cover 30,000 deg^, with two visits per year 
with rUg < 23.54, over 3 yr, resulting in a detection rate 
of ~130 yr~^. The results of the simulation were fitted 
with a polynomial for clarity. The LSST 20,000 deg^ 
survey will increase these rates by 2 orders of magnitude, 
since it has a similar depth and cadence to the MDS, with 
individual visits with nig < 24.5, but the factor of 400 
increase in survey area that will yield 6000 events yr^^! 

The next generation of wide-field optical synoptic sur- 
veys will increase the number of tidal disruption flare de- 
tections by orders of magnitude, allowing for the study 
of the ensemble properties and rates of tidal disruption 
events as a function of galaxy type. The rates of stellar 
disruption events are sensitive to the detailed structure 
of the galaxy nuclei, and can be used to test dynamical 
models for galaxies that are too distant for detailed spa- 
tially resolved observations. The light curves of individ- 
ual tidal disruption events can be used as a direct probe 
of the central black hole masses, and a means to map 
the dormant black hole mass function in normal galaxies. 
The SEDs of events may also be able to probe the spin 
of black holes, since the characteristic radius of the emis- 
sion can be compared to the smallest stable particle orbit, 
which is directly dependent on the spin. Although the de- 
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Fig. 27. — Number of events predicted to be detected in the g 
band per year within a volume out to redshift z = 0.5 for fc = 3 as 
a function of black hole mass in the Pan-STARRS MDS and the 
Stt Survey. 

tection sensitivities of these surveys are "tuned" to black 
hole masses of a limited mass range, (1 — 3) x 10 ''M©, 
they will be a sensitive measure for deviations from the 
locally established scaling relations between the mass of 
the black hole and their host galaxy properties. Evidence 
for such deviations has already been seen for a sample of 
14 Se yfert gala xies at z - 0.37 with A/bh = 10^ - lO^M© 
bv Woo et al. I (^2006). Large samples of tidal disruption 
event detections will yield a fruitful data set for testing 
models for black hole growth and galaxy evolution. 
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TABLE 1 
GALEX Yearly Coadds 



Field 


a (J2000) 


<5 (J2000) 


E(B-V) 


Band 




Exposure Time (ks) 














2003 


2004 


2005 


2006 


2007 


XMM/LSS 


36.3610 


-4.4800 


0.027 


NUV 


17.8 


1.1 


27.5 


20.0 












FUV 


17.8 


1.1 


22.1 


20.0 




C09M0S 


149.8690 


2.4558 


0.018 


NUV 




47.5 




8.6 


20.0 










FUV 




41.5 




8.6 


20.0 


GROTH 


214.9918 


52.7817 


0.010 


NUV 


16.1 


1.5 


102.4 


60.1 


15.7 










FUV 


16.1 


1.5 


40.7 


42.7 


15.7 


CPHTLS D4 


333.7830 


-17.9350 


0.027 


NUV 


24.3 


5.0 


27.6 


18.2 












FUV 


24.3 


5.0 


5.6 


18.1 





TABLE 2 

X'AiiiADLE Sources with Optical Matctieh ix XMM/LSS (D1) 



o 



No. 




a(JzUOUj 




2003 


2004 


r U V 
2005 


ZUOO 2007 


a 






2003 


2004 


l\U V 
2005 


2006 2007 


Dl 


1 


J022409. 7-044007 


36.0406 


-4.6687 


23.72±0.13 


>22.8 


22.61±0.05 


22.93±0.06 


8 


8 


F 


22 


72±0.05 


>22.7 


22.01±0.03 


22 


22±0.03 


Dl 


2 


J022409. 7-043405 


36.0406 


-4.5682 


22.61±0.06 


23.00±0.31 


24.00±0.13 


23.18±0.08 


9 


4 




21 


53±0.02 


22,24±0.15 


22.44±0.04 


21 


98±0.03 


Dl 


3 


J022424. 2-043230 


36.1009 


-4.5418 


22.44±0.05 


23.22±0.35 


23.52±0.09 


23.23±0.08 


8 


1 




20 


03±0.01 


20.89±0.06 


20.98±0.01 


21 


OOiO.Ol 


Dl 


4 


J022434. 7-045350 


36.1444 


-4.8973 


>24.3 


>22.8 


23.91±0.12 


>24.3 


6 





T 




>24.2 


>22,7 


20.38±0.01 




>24.3 


Dl 


5 


J022449. 2-041802 


36.2052 


-4.3006 


22.79±0.07 


>22.8 


24.04±0.13 


24.80±0.24 


16 


1 


F 


22 


12±0.03 


>22,7 


22.75d=0.05 


23 


16±0.07 


Dl 


6 


J022449. 9-043025 


36.2078 


-4.5072 


22.92±0.07 


22.34d=0.20 


23.54±0.09 


23.89±0.12 


6 


3 




21 


42±0.02 


20.93±0.06 


21.71±0.02 


21 


93±0.03 


Dl 


7 


J022451. 2-040527 


36.2132 


-4.0911 


22.47±0.06 


23.16±0.35 


22.55±0.05 


23.16±0.08 


5 


2 




21 


27±0.02 


21.50±0.09 


21.34±0.02 


21 


78±0.02 


Dl 


8 


J022503. 2-040539 


36.2635 


-4.0944 


22.77±0.07 


23.35±0.39 


21.43±0.03 


21.40±0.03 


14 


7 




21 


46±0.02 


22.20±0.14 


20.65±0.01 


20 


58±0.01 


Dl 


9 


,1022517.0-043258 


36.3207 


-4.5497 


>24.3 


23.23±0.36 


23.41±0.08 


24.41±0.18 


7 


7 


T 




>24.2 


23.05±0.27 


23.32±0.07 


24 


00±0.13 


Dl 


10 


J022518. 1-043155 


36.3254 


-4.5320 


23.86±0.14 


>22.8 


>24.4 


>24.3 


6 


9 


F 


23 


47±0.10 


>22.7 


>24.4 




>24.3 


Dl- 


11 


J022526. 5-040305 


36.3604 


-4.0515 


>24.3 


>22.8 


24.38±0.17 


23.65±0.11 


5 


5 


F 




>24.2 


>22.7 


22.79±0.05 


22 


06±0.03 


Dl- 


12 


J022527. 5-035954 


36.3647 


-3.9984 


23.89±0.15 


23.17±0.35 


23.64±0.10 


24.72±0.23 


5 


7 




22 


69±0.05 


22.09±0.13 


22.82±0.05 


23 


75±0.11 


Dl- 


13 


J022539.4-042228 


36.4141 


-4.3746 


22.12±0.04 


23.23±0.36 


22.45±0.05 


24.55±0.20 


25 







21 


25±0.02 


21.97±0.12 


21.76±0.02 


22 


86±0.05 


Dl- 


14 


J022544.4-040102 


36.4350 


-4.0173 


22.65±0.06 


23.28±0.38 


23.18±0.07 


23.93±0.13 


8 


6 


F 


21 


85±0.03 


21.81±0.11 


22.12±0.03 


22 


54±0.04 


Dl- 


15 


J022552.2-042249 


36.4674 


-4.3803 


22.56±0.06 


22.12±0.18 


22.76±0.05 


23.55±0.10 


7 


4 




21 


31±0.02 


21.05±0.07 


21.50±0.02 


21 


91±0.03 


Dl 


16 


J022555.4-043918 


36.4810 


-4.6551 


21.96±0.04 


22.36±0.21 
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20 
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20.63±0.05 


20.56±0.01 


21 


09±0.02 


^ FUV flares that fade monotonically, have a steady flux before the flare, and do not fade below the steady-state flux are labeled with an 


'F'. Transient FUV flares that fade monotonically, and are not detected before the start of 



the flare are labeled with a 'T'. 



TABLE 3 

UV Variable SorBCES with Optkwl Matciies in COSMOS (D2) 



INO. 


GALEX ID 
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^ FUV flares that fade monotonically, have a steady flux before the flare, and do not fade below the steady-state flux are labeled 


with 


an 


'F'. 


Transient FUV flares that fade monotonically, and 


are not detected 



before the start of the flare are labeled with a 'T'. 
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.yoitu. Uo 


01 

zl 


Qfi-Un no 
.oDitU. UZ 


OO 
ZZ 


Q 1 -Ln Of; 
.vSlitU.UO 


r: 
0. 


D 




on 

zU. 


QQ_LO 01 
.OOltU.Ul 


oo /IT-LO O/l 

zU.4 / itU.U4 


on fi/1 -LO 01 
zU. D4itU.Ul 


on 70-L0 OI 

ZU. ( zitU.Ul 


OI 

zl. 


1 fi-LO OO 

.iDitU.Uz 


RQ 


1 o 

-ly 


11/1 ori/i 


9_L f^Qflfil A 

.Z-j-OoUolD 


91 

ZlO. 


9ni n 

.ZUlU 


Oo. 


.lo / y 




\9/1 
>Z4.0 




\OQ o 


99 
ZZ. 


fiO-UO 0/1 

.oUztU.U4 


9Q 

Zo 


fifi-UO 07 




\0/1 f; 
^Z4.0 


1 Q 
lo. 


1 

. 1 


1 




\0/1 A 
^Z4.4 


\OQ 1 
>Zo. 1 


09 /10-L0 OQ 
ZZ. 4ZltU.Uo 


OQ 07-L0 Of; 
Zo.Z ( itU.UO 




\0/1 A 
^Z4.4 


D3- 


■20 


J142056, 


,9-f524830 


215, 


,2370 


52, 


,8084 


21 


.83+0.04 


22 


.33ib0.17 


22 


.33+0.03 


21 


.66+0.02 


22 


.34+0.05 


6, 


.2 




20. 


,78+0.01 


21.00+0.05 


21.08+0.01 


20.52+0.01 


21, 


,07+0.02 


D3- 


■21 


J142107, 


,7+530319 


215, 


.2821 


53, 


,0553 


23 


.81+0.15 




>23.2 


22 


.75+0.04 


23 


.33+0.06 


24 


.09+0.17 


8, 


.0 


F 


22. 


.40+0.04 


>23.1 


22.02+0.02 


22.22+0.02 




,70+0.05 


D3- 


■22 


J142112, 


,3+524147 


215, 


,3012 


52, 


,6965 


21 


.01+0.03 


21 


.61it0.12 


21. 


.67+0.02 


21 


.85+0.02 


21 


.61+0.03 


6, 


.0 




20. 


.30+0.01 


20.70+0.05 


20.76+0.01 


20.77+0.01 


m 


,64+0.01 


D3- 


■23 


J142113, 


,9+521747 


215, 


.3081 


52, 


,2966 


22 


.50+0.06 


22 


.32ib0.17 


22 


.34+0.03 


22 


.06+0.03 


23 


.06+0.08 


7, 


.0 




22. 


.36+0.04 


22.23+0.12 


22.23+0.02 


21.96+0.02 


^2, 


,90+0.06 


D3- 


■24 


J142115, 


,7+523156 


215, 


.3154 


52, 


,5323 


23 


.89+0.16 




>23.2 


23 


.86+0.09 


23 


.14+0.05 


24 


.01+0.16 


5, 


.5 




22. 


.47+0.05 


>23.1 


22.68+0.03 


22.19+0.02 




,96+0.06 


D3- 


■25 


J142125, 


,0+521800 


215, 


.3542 


52, 


,3001 




>24.5 




>23.2 


25 


.49+0.31 


22 


.29+0.03 


25 


.09+0.37 


33, 


.6 


F 




>24.4 


>23.1 


23.04+0.04 


21.24+0.01 


^3, 


,13+0.07 


D3- 


■26 


J142133, 


,9+530245 


215, 


.3912 


53, 


,0460 


23 


.86+0.16 


23 


.51ib0.38 


23 


.35+0.06 


23 


.23+0.05 


22 


.53+0.06 


5, 


.9 




22 


.39+0.04 


22.00+0.10 


22.07+0.02 


22.04+0.02 


esi 


,61+0.02 


D3- 


■27 


J142134, 


,2+523429 


215, 


.3925 


52, 


,5749 


23. 


.75+0.14 




>23.2 


25 


.10+0.23 


23 


.91+0.08 


23 


.75+0.13 


5, 


.5 




21. 


.15+0.02 


>23.1 


21.67+0.02 


21.21+0.01 


^1, 


.15+0.02 


D3- 


■28 


J142135, 


,9+523139 


215, 


.3996 


52, 


,5277 


20 


.49+0.02 


21 


.09±0.09 


20 


.07+0.01 


20 


.56+0.01 


21 


.03+0.03 


7, 


.5 




20. 


.28+0.01 


20.77+0.05 


19.95+0.01 


20.17+0.01 


20, 


.63+0.01 


D3- 


■29 


J142151. 


,1+524951 


215, 


.4629 


52, 


,8309 


22 


.24+0.05 




>23.2 


21 


.67+0.02 


22 


.19+0.03 


22 


.21+0.05 


5, 


.3 




21. 


.85+0.03 


>23.1 


21.44+0.01 


21.84+0.02 


21, 


.84+0.03 


D3- 


■30 


J142206. 


,8+524957 


215, 


.5285 


52, 


,8327 


23 


.23+0.10 




>23.2 




>25.0 


23 


.66+0.07 


24, 


.16+0.17 


7, 


.3 




22. 


.68+0.05 


>23.1 


>25.4 


22.79+0.03 


23, 


.11+0.07 


D3- 


•31 


J142209, 


,1+530559 


215, 


.5381 


53, 


,1000 


20 


.94+0.02 


21 


.llibO.09 


20 


.88+0.02 


21 


.67+0.02 


21 


.38+0.03 


6, 


.4 




19. 


.81+0.01 


19.93+0.03 


19.85+0.01 


20.26+0.01 


20, 


.15+0.01 


D3- 


■32 


J142213 


,5+530520 


215 


,6811 


53, 


,0890 


22 


,.39+0.06 




>23.2 


22 


.09+0.03 


21 


.65+0.02 


23 


.21+0.09 


13 


,7 




21, 


,40+0.02 


->23.1 


21.06+0.01 


20.65+0.01 


21 


,72+0.03 



^ PUV flares that fade monotonically, have a steady flux before the flare, and do not fade below the steady-state flux are labeled with an 'F'. Transient FUV flares that fade monotonically, and are not detected before the start of the flare are 
labeled with a 'T'. 



TABLE 5 

UV Variable Sources with Optical Matches in CFHTLS D4 



No. 


GALEX ID 


a(J2000) 


5(J2000) 






FUV 










NUV 














2003 


2004 


2005 


2006 2007 




2003 


2004 


2005 


2006 2007 


D4-1 


J221336.1- 


180531 


333.4003 


-18.0920 


21.25±0.02 


21.15±0.05 


22.28±0.08 


21.34±0.03 


12.7 


20.20±0.01 


20.07±0.02 


20.97±0.03 


20.41±0.01 


D4-2 


J221409.7- 


175652 


333.5403 


-17.9480 


21.25±0.02 


21.15±0.05 


21.06±0.04 


20.72±0.02 


5.4 


20.31±0.01 


20.31±0.02 


20.13±0.02 


19.94±0.01 


D4-3 


J221413.0- 


172909 


333.5542 


-17.4860 


22.66±0.05 


22.41±0.10 


21.65±0.06 


23.16±0.08 


8.7 


21.18±0.02 


20.90±0.03 


20.33±0.02 


21.69±0.02 


D4-4 


J221429.5- 


173007 


333.6230 


-17.5020 


19.29±0.01 


19.28±0.02 


19.71±0.02 


19.66±0.01 


5.4 


18.37±0.00 


18.35±0.01 


18.65±0.01 


18.65±0.00 


D4-5 


J221432.7- 


180158 


333.6363 


-18.0.330 


21.41±0.03 


21.80±0.07 


21.11±0.04 


22.01±0.04 


7.2 


20.19±0.01 


20.51±0.02 


19.84±0.01 


20.40±0.01 


D4-6 


J221433.8- 


173007 


333.6408 


-17.5020 


22.91±0.06 


23.55it0.20 


23.26±0.15 


22.57±0.05 


5.5 


22.28it0.03 


22.49it0.08 


22.38±0.07 


22.00±0.03 


D4-7 


J221508.9- 


175602 


333.7869 


-17.9340 


23.43±0.08 


>23.6 


>23.7 


>24.3 


6.5 F 


22.50±0.04 


>23.5 


>24.4 


>24.2 


D4-8 


J221639.6- 


174427 


334.1651 


-17.7410 


22.83±0.06 


21.28±0.05 


21.88±0.07 


21.54±0.03 


12.2 


21.37±0.02 


20.41±0.02 


20.96±0.03 


20.56±0.01 


D4-9 


.1221617.8- 


1751.51 


331.1991 


-17.91.o0 


22.11±0,01 


21.91±0.07 


22.35±0.09 


22.82±0.06 


'',.2 


21.23±0,02 


20.90±0,03 


21.19±0,03 


21..o3±0,02 



^ FUV flares that fade monotonically, have a steady flux before the flare, and do not fade below the steady-state flux are labeled with an 'F'. Transient FUV flares that fade monotonically, and are not detected before the start 
of the flare are labeled with a 'T'. 



24 Gezari et al. 



TABLE 6 

Optical Matches to UV Variable Sources in XMM/LSS (D1) 



No. 


GALEX ID 




u 






9 




r 




i 






z 




^1/2 

(arcsec) 


Sep 
(arcsec) 


flag^ 


class'' z 


Dl- 


■1 


J022409. 


.7-044007 


22, 


.lllibO. 


.007 


21, 


.288±0. 


.003 


20.286±0.002 


19. 


.665±0, 


.001 


19, 


.298±0. 


.002 


1.043 


0.6 


o 


GAL 


Dl- 


■2 


J022409. 


.7-043405 


21, 


.982±0. 


.003 


21, 


.703±0. 


.002 


21.359±0.002 


21. 


.298±0, 


.002 


21, 


.348±0. 


.008 


0.461 


1.0 


X 


QSO 


Dl- 


■3 


J022424, 


.2-043230 


19 


.550±0. 


.000 


19, 


.192±0. 


.000 


19.171±0.000 


18. 


.923±0, 


.000 


18, 


.972±0. 


.001 


0.452 


1.0 


X 


QSO 


Dl- 


■4 


J022434. 


.7-045350 


18, 


.llliO. 


.000 


17, 


.102±0. 


.000 


16.703±0.000 


16. 


.612±0, 


.000 


16, 


.420±0. 


.000 


0.500 


2.2 




STAR 


Dl- 


■5 


J022449. 


.2-041802 


22, 


.476±0. 


.012 


21, 


.750±0. 


.005 


21.032±0.003 


20. 


.419±0, 


.002 


20, 


.132±0. 


.006 


0.703 


0.2 


X 


AGN 


Dl- 


■6 


J022449, 


.9-043025 


20 


.613±0. 


.001 


20 


.319±0, 


.001 


20.002±0.001 


19. 


.982±0 


.001 


19 


.850±0, 


.002 


0.472 


0.2 


X 


QSO 


Dl- 


-7 


J022451, 


,2-040527 


20, 


.754±0, 


,001 


20, 


.521±0, 


,001 


20.491±0.001 


20, 


,276±0, 


.001 


20, 


.278±0, 


,004 


0.458 


0.6 




QSO 


Dl- 


-8 


J022503, 


,2-040539 


21 


.837±0, 


,004 


21 


.090±0, 


,002 


20.980±0.002 


20, 


,791±0, 


.002 


20, 


.458±0, 


,005 


0.564 


0.9 


X 


QSO 


Dl- 


■9 


J022517, 


,0-043258 


23 


.243±0, 


,015 


22 


.106±0, 


,005 


20.950±0.002 


20, 


,369±0 


.001 


20 


.064±0, 


,004 


0.636 


0.1 


o 


GAL 


Dl- 


■10 


,1022518, 


,1-043155 


24, 


.450±0, 


,039 


24, 


.272±0, 


,032 


23.360± 0.016 


22, 


,917±0 


.013 


23 


.059±0, 


,051 


0.712 


0.6 


o 


GAL 


Dl- 


■11 


J022526, 


,5-040305 


22 


.736±0, 


,007 


22 


.636±0, 


,006 


22.114± 0.004 


21. 


,947±0 


.004 


21 


.362±0, 


,009 


0.528 


0.2 


o 


QSO 


Dl- 


■12 


J022527, 


,5-035954 


21, 


.528±0, 


,003 


21, 


.254±0, 


,002 


21.225± 0.003 


20, 


,951±0, 


.002 


20, 


.852±0, 


,007 


0.439 


0.3 




QSO 


Dl- 


■13 


J022539, 


,4-042228 


21 


.955±0, 


,004 


21 


.457±0, 


,002 


21.269± U.002 


20. 


,953±0 


.002 


20 


.567±0, 


,005 


0.508 


0.7 




QSO 


Dl- 


■14 


J022544, 


,4-040102 


22 


.042±0, 


,005 


21 


.637±0, 


,004 


21.396± 0.003 


20, 


,973±0, 


.003 


20, 


.700±0, 


,007 


0.623 


0.4 




GAL 


Dl- 


■15 


J022552, 


,2-042249 


21, 


,417±0, 


,002 


21, 


.056±0, 


,001 


20.982± 0.002 


21, 


,015±0, 


.002 


20, 


.697±0, 


,005 


0.437 


0.5 


X 


QSO 


Dl- 


■16 


J022555, 


,4-043918 


21, 


,069±0, 


,002 


20, 


.823±0, 


,001 


20.568± 0.001 


20, 


,517±0, 


.001 


20, 


.273±0, 


,003 


0.437 


0.3 


X 


QSO 


Dl- 


■17 


J022556, 


,9-045853 


21, 


.484±0, 


,003 


21, 


.109±0, 


,002 


20.790± 0.002 


20, 


,751±0, 


.002 


20, 


.591±0, 


.006 


0.515 


0.7 




QSO 


Dl- 


■18 


J022558, 


.1-045720 


20, 


.962±0. 


.002 


20, 


.677±0. 


.001 


20.386± 0.001 


20. 


.378±0, 


.001 


20, 


.238±0. 


.004 


0.497 


0.6 




QSO 


Dl- 


■19 


J022606. 


.7-045723 


21, 


.004±0. 


.002 


20, 


.709±0. 


.001 


20.629± 0.001 


20. 


.264±0, 


.001 


20, 


.026±0. 


.003 


0.560 


0.2 


ox 


QSO 


Dl- 


■20 


J022639. 


.8-042004 


21, 


.525±0. 


.002 


21, 


.640d=0. 


.003 


21.430± 0.003 


21. 


.130±0, 


.002 


21, 


.111±0. 


.008 


0.419 


0.2 


X 


QSO 


Dl- 


■21 


J022647. 


.8-043135 


22, 


.720±0. 


.007 


22, 


.518d=0. 


.006 


22.052± 0.004 


21. 


.851±0, 


.004 


21, 


.410±0. 


.010 


0.502 


0.7 




QSO 


Dl- 


■22 


J022701, 


.5-040913 


21, 


.389±0. 


.003 


21, 


.021±0. 


.002 


20.701± 0.002 


20. 


.227±0, 


.001 


19 


.964±0. 


.003 


0.731 


0.5 


OX 


AGN 


Dl- 


■23 


J022716, 


.1-044539 


20, 


.939±0. 


.002 


20, 


.320±0, 


.001 


20.309± 0.001 


20. 


.152±0, 


.001 


19, 


.856±0, 


.003 


0.543 


0.5 


X 


QSO 


a „ 

Sources 


that are optically variable are flagged 


with 


an 'o 


', sources 


i detected as a hard X-ray poiin 


; source 


are flagged with an 


'x', sources with 


an optical AGN spectrum are flagged 



with an 's', and sources with a UV quasar spectrum are flagged with a 'g'. 

^ Optically unresolved sources (rj^y2 0^6) with g — r < 0.6 and u — g < 1 are classified as quasars (QSO). Optically unresolved sources with u — g > 1.75 or g — r < 0.6 and u — g > 1 
are classified as stars (STAR). Optically resolved sources (^1/2 0^'6) are classified as galaxies (GAL), and optically resolved sources with a hard X-ray detection or an active galaxy 
spectrum are classified as active galactic nuclei (AGN). 



TABLE 7 

Optkwl Matches to UV Variable Soprceh ix COSMOS (D'2) 



No. 


AT TFi 




U 






9 




r 




2 






z 




^1/2 

(arcsec) 


Sep 
(arcsec) 


nag 


1 b 

class 


Z 


D2- 




J095829. 3+021542 


21 


.665±0, 


Oil 


21. 


.297+0. 


.005 


21.000+0.005 


20. 


.617+0. 


.004 


20 


.284+0, 


.008 


0.586 


0.7 




QSO 






D2- 


.2 


J095848. 9+023441 


20, 


.421+0, 


.002 


20. 


.849+0. 


,002 


20.606+0.002 


20, 


,344+0. 


,002 


20 


.456+0, 


.005 


0.498 


0.5 


xs 


QSO 


I 


551 


D2- 


.3 


J095858. 6+021459 


19, 


.486+0, 


.006 


18. 


.054+0. 


,001 


17.253+0.001 


16, 


.779+0, 


,000 


16 


.528+0, 


.001 


2.052 


0.4 


xsg 


AGN 


0. 


132 


D2- 


.4 


J095921. 5+024029 


19, 


.504+0, 


.003 


18, 


.614+0. 


,001 


17.794+0.001 


17, 


.308+0. 


,000 


17, 


.002+0, 


.001 


1.949 


2.0 


SET 


AGN 


0. 


.260 


D2- 


.5 


J095947.0+022209 


21 


.304+0, 


.004 


20, 


.749+0. 


,002 


20.697+0.002 


20, 


.763+0. 


,002 


20, 


.554+0, 


.006 


0.472 


0.7 


xs 


QSO 


0. 


909 


D2- 


-6 


J095954. 8+021706 


17, 


.345+0, 


.001 


15, 


.458+0. 


,000 


14.634+0.000 


14, 


.320+0. 


,000 


13, 


.625+0, 


.000 


0.500 


1.3 




STAR 






D2- 


.7 


J095958. 5+021530 


21 


.719+0, 


.013 


21, 


.501+0. 


,004 


21.043+0.003 


20, 


.145+0. 


,002 


19, 


.794+0, 


.003 


0.731 


0.8 


X 


AGN 






D2- 


-8 


J095958. 9+020023 


23, 


.167+0, 


.020 


22, 


.621+0. 


,007 


22.146+0.006 


21, 


.831+0. 


,005 


21 


.459+0, 


.010 


0.515 


1.1 





QSO 






D2- 


-9 


J100002. 0+024216 


22, 


.575+0, 


.018 


21, 


.880+0. 


,006 


21.410+0.004 


20, 


.771+0. 


,003 


20, 


.464+0, 


.007 


0.655 


0.4 




GAL 






D2- 


■10 


J100012. 9+023522 


19, 


.454+0, 


.001 


18, 


.941+0. 


,000 


18.868+ 0.000 


18, 


.605+0. 


,000 


18, 


.473+0, 


.001 


0.500 


0.4 




QSO 


0. 


.699 


D2- 


-11 


J100014. 9+022717 


21 


.882+0, 


.011 


21, 


.114+0. 


,004 


20.964+ 0.005 


20, 


.638+0. 


,004 


20, 


.356+0, 


.007 


0.658 


0.3 


GAL 






D2- 


■12 


J100017.6+020011 


22, 


.148+0, 


.010 


21, 


.993+0. 


,005 


21.095+ 0.003 


20, 


.664+0. 


.002 


20, 


.299+0, 


.005 


0.569 


1.2 










D2- 


■13 


J100024.7+023149 


18, 


.915+0, 


.000 


18, 


.814+0. 


.000 


18.614+ 0.000 


18, 


.663+0. 


.000 


18, 


.748+0, 


.001 


0.456 


0.5 


xsg 


QSO 


1. 


.318 


D2- 


■14 


J100029. 7+022129 


21 


.754+0, 


.007 


21, 


.023+0. 


.002 


20.775+ 0.002 


20. 


.217+0. 


.002 


19, 


.923+0, 


.004 


0.729 


0.2 





GAL 






D2- 


■15 


J100035.0+020235 


21, 


.079+0. 


.004 


20. 


.862+0. 


.002 


20.549+ 0.002 


20. 


.478+0. 


.002 


20, 


.402+0. 


.005 


0.476 


0.8 


xs 


QSO 


1. 


.177 


D2- 


■16 


J100046.8+020404 


20, 


.701+0. 


.004 


20. 


.386+0. 


.002 


20.043+ 0.002 


19. 


.549+0. 


.001 


19, 


.324+0. 


.003 


0.831 


0.7 


oxs 


AGN 


0. 


.554 


D2- 


■17 


J100055.4+023441 


20, 


.802+0. 


.007 


20. 


.558+0. 


.001 


20.187+ 0.001 


20. 


.143+0. 


.001 


20, 


.135+0. 


.004 


0.458 


0.5 


xs 


QSO 


1. 


.403 


D2- 


■18 


J100104.2+023349 


25, 


.030+0. 


.273 


24. 


.565+0. 


.041 


24.468+ 0.047 


24. 


.196+0. 


.043 


23 


.765+0. 


.079 


0.497 


1.7 




QSO 






D2- 


■19 


J100104.3+023402 


22, 


.452+0, 


.046 


21, 


.838+0. 


,006 


21.150+ 0.004 


20, 


.882+0. 


,004 


20, 


.739+0, 


.009 


0.900 


0.8 




GAL 






D2- 


■20 


J100104.8+023329 


25, 


.831+0, 


.566 


25, 


.902+0. 


141 


25.842+ 0.170 


25, 


.798+0. 


,189 


24, 


.989+0, 


.246 


0.374 


2.9 




QSO 






D2- 


■21 


J100105.4+023323 


22, 


.652+0, 


.051 


21, 


.383+0. 


,004 


20.512+ 0.002 


20, 


.082+0. 


,002 


19, 


.800+0, 


.004 


0.683 


1.0 




ART 






D2- 


■22 


J100116.3+023607 


21, 


.666+0, 


.005 


20, 


.920+0. 


,002 


20.785+ 0.002 


20, 


.608+0, 


.002 


20, 


.308+0, 


.004 


0.582 


0.4 




QSO 






D2- 


•23 


J100124.0+021445 


22, 


.887+0, 


.014 


22. 


.241+0. 


.005 


22.120+ 0.005 


21. 


.868+0. 


.005 


21, 


.524+0, 


.011 


0.484 


0.3 


xs 


QSO 


0. 


.894 


^ Sources 


that are optically variable 


are : 


Flagged w 


'ith an 


'0', 


sources detected 


, as a hard X-ray point source are 


i flagged with an 'x' 


, sources with an 


Optical AGN spectrum 


are flagged with an 



's', and sources with a UV quasar spectrum are flagged with a 'g'. 

^ Optically unresolved sources (X\i'2, < 0^6) with g — r < 0.6 and u — g < 1 are classified as quasars (QSO). Optically unresolved sources with tx — g > 1.75 or g — r < 0.6 and tx — g > 1 are 
classified as stars (STAR). Optically resolved sources (rj^^2 ^ 0^6) are classified as galaxies (GAL), and optically resolved sources with a hard X-ray detection or an active galaxy spectrum 
are classified as active galactic nuclei (AGN). 



TABLE 8 

Optical Matches to \ 'V Nariahle Socrceh ix (4H()TH (D?>) 



IMO. 


AT TPi 




U 






9 




r 






2 






z 




^1/2 


Sep 


flar 




1 b 

class 


Z 


































(arcsec) 


(arcsec) 












D3-1 


J141714.4+525130 


26 


.621+0. 


.332 


25. 


.568+0, 


.071 


25.243+0.060 


24. 


.982+0 


.060 


26 


.445+0. 


.787 


0.439 


1.7 






STAR 






D3-2 


J141715. 2+530303 


20, 


.979+0. 


.002 


20. 


.915+0, 


.001 


20.865+0.001 


20. 


.784+0, 


.002 


20, 


.550+0. 


.004 


0.472 


0.6 


o 




QSO 






D3-3 


J141721. 2+530916 


22, 


.108+0. 


.013 


21. 


.333+0, 


.004 


20.061+0.001 


19. 


.503+0, 


.001 


19, 


.162+0. 


.002 


0.945 


0.4 


o 




GAL 






D3-4 


J141724. 6+523025 


20, 


.033+0. 


.001 


19. 


.686+0, 


.001 


19.586+0.001 


19. 


.324+0, 


.001 


19, 


.124+0. 


.001 


0.500 


0.7 


X 


6 


QSO 


0, 


.482 


D3-5 


J141734. 9+522811 


21 


.988+0. 


.006 


21. 


.951+0, 


.003 


21.410+0.002 


21. 


.338+0, 


.003 


20, 


.841+0. 


.006 


0.508 


0.8 


X 




QSO 






D3-6 


J141801. 5+525201 


20, 


.261+0. 


.001 


20. 


.162+0, 


.001 


19.991+0.001 


19. 


.963+0, 


.001 


19, 


.849+0. 


.003 


0.474 


0.3 


6 




QSO 


1, 


.159 


D3-7 


J141830. 3+522213 


21 


.520+0. 


.004 


20. 


.740+0, 


.001 


20.713+0.001 


20. 


.546+0, 


.001 


20, 


.243+0. 


.004 


0.523 


0.6 


X 


& 


QSO 





.825 


D3-8 


J141833. 5+530733 


22, 


.957+0. 


.016 


22. 


.365+0, 


.006 


21.726+0.004 


20. 


.969+0, 


.002 


20, 


.666+0. 


.006 


0.642 


0.3 


o 


GAL 






D3-9 


J141834. 9+524205 


25, 


.825+0. 


.196 


25. 


.581+0, 


.075 


24.649+0.035 


24. 


.319+0, 


.033 


23, 


.860+0. 


.079 


0.458 


1.6 












D3-10 


J141903. 9+530856 


22, 


.533+0. 


.009 


21. 


.925+0, 


.003 


21.692+ 


0.003 


21. 


.411+0, 


.003 


21 


.484+0. 


.011 


0.441 


0.4 






QSO 






D3-11 


J141905. 2+522528 


19, 


.560+0. 


.001 


19. 


.539+0, 


.001 


19.503+ 


0.001 


19. 


.322+0, 


.001 


19, 


.313+0. 


.002 


0.437 


0.7 


Sfi 




QSO 


1, 


.606 


D3-12 


J141918. 1+524158 


20, 


.426+0. 


.001 


20. 


.322+0, 


.001 


20.297+ 


0.001 


20. 


.149+0, 


.001 


19, 


.930+0. 


.003 


0.430 


0.3 


g 




QSO 


2, 


.034 


D3-13 


J141929. 8+525206 


22, 


.741+0. 


.025 


21. 


.363+0, 


.004 


20.108+ 


0.002 


19. 


.634+0, 


.001 


19, 


.349+0. 


.004 


0.625 


0.4 


o 




GAL 






D3-14 


J141936. 1+521359 


23, 


.179+0. 


.021 


22. 


.141+0, 


.005 


21.374+ 


0.003 


20. 


.985+0, 


.002 


20, 


.762+0. 


.007 


0.688 


1.5 






ERR 






D3-15 


J141946. 4+525943 


17. 


.619+0. 


.000 


15. 


.695+0. 


.000 


14.828+ 


0.000 


14. 


.138+0. 


.000 


13 


.443+0. 


.000 


0.500 


0.5 






STAR 






D3-16 


J141951. 5+524205 


22. 


.711+0. 


.016 


22. 


.066+0. 


.006 


21.404+ 


0.003 


20. 


.678+0. 


.002 


20. 


.405+0. 


.006 


0.766 


0.6 


o 




GAL 






D3-17 


J141959. 8+524245 


22. 


.253+0. 


.008 


21. 


.403+0. 


.002 


20.649+ 


0.001 


20. 


.255+0. 


.001 


20. 


.039+0. 


.003 


0.508 


0.2 


o 










D3-18 


J142039. 5+521929 


19 


.519+0. 


.001 


19. 


.323+0 


.000 


18.985+ 


0.000 


18. 


.665+0, 


.000 


18. 


.587+0. 


.001 


0.437 


0.6 




> 


QSO 


1. 


.590 


D3-19 


J142048. 2+530816 


23, 


.139+0. 


.013 


23. 


.791+0, 


.014 


23.672+ 


0.015 


23. 


.489+0, 


.015 


22, 


.899+0. 


.030 


0.450 


0.3 






QSO 






D3-20 


J142056. 9+524830 


20, 


.653+0. 


.002 


20. 


.548+0, 


.001 


20.187+ 


0.001 


20. 


.127+0, 


.001 


20, 


.074+0. 


.004 


0.458 


0.3 


& 




QSO 


1. 


.199 


D3-21 


J142107.7+530319 


22, 


.383+0. 


.009 


21. 


.871+0, 


.004 


21.529+ 


0.003 


20. 


.958+0, 


.002 


20, 


.639+0. 


.005 


0.629 


0.9 


O 




GAL 






D3-22 


J142112.3+524147 


21, 


.111+0. 


.003 


20. 


.446+0, 


.001 


20.369+ 


0.001 


20. 


.298+0 


.001 


20, 


.085+0, 


.003 


0.506 


0.1 


g 




QSO 


0, 


.847 


D3-23 


J142113.9+521747 


19, 


.934+0, 


.001 


19. 


.715+0 


.001 


19.703+ 


0.001 


19. 


.565+0 


.001 


19, 


.451+0, 


.002 


0.459 


0.6 




QSO 






D3-24 


J142115.7+523156 


21, 


.674+0, 


.004 


21. 


.301+0 


.002 


21.062+ 


0.002 


20. 


.933+0 


.002 


20, 


.913+0, 


.007 


0.428 


0.8 






QSO 






D3-25 


J142125.0+521800 


22, 


.971+0, 


.013 


22. 


.523+0 


.005 


22.299+ 


0.005 


21. 


.989+0 


.005 


21, 


.483+0, 


.010 


0.497 


0.5 






QSO 






D3-26 


J142133.9+530245 


22, 


.141+0, 


.006 


21. 


.702+0 


.003 


21.556+ 


0.003 


21. 


.350+0 


.003 


21, 


.040+0, 


.007 


0.487 


0.2 






QSO 






D3-27 


J142134.2+523429 


21, 


.194+0, 


.003 


21. 


.007+0 


.001 


20.688+ 


0.001 


20. 


.611+0 


.002 


20, 


.435+0, 


.004 


0.465 


0.5 


g 




QSO 


1, 


.199 


D3-28 


J142135.9+523139 


19, 


.601+0, 


.002 


19. 


.050+0 


.001 


18.298+ 


0.001 


17. 


.859+0 


.000 


17, 


.624+0, 


.001 


1.326 


0.6 


o 


sg 


AGN 


1, 


.249 


D3-29 


J142151. 1+524951 


20 


.786+0, 


.005 


19. 


.615+0 


.001 


18.818+ 


0.001 


18. 


.355+0 


.001 


18 


.111+0, 


.001 


1.196 


0.3 


o 




GAL 






D3-30 


J142206. 8+524957 


22 


.740+0, 


.009 


22. 


.076+0 


.003 


22.139+ 


0.004 


21. 


.690+0 


.003 


21 


.853+0, 


.014 


0.439 


0.8 






QSO 






D3-31 


J142209.1+530559 


19, 


.907+0. 


.001 


19. 


.515+0, 


.001 


19.429+ 


0.001 


19. 


.207+0 


.001 


18, 


.902+0. 


.001 


0.597 


0.1 


g 


QSO 


0, 


.760 


D3-32 


J142243.5+530520 


20 


.401+0. 


.005 


20. 


.407+0, 


.003 


20.122+ 


0.003 


20. 


.053+0 


.003 


19 


.938+0. 


.011 


0.621 


0.3 






GAL 






^ Sources 


tliat arc optically variable 




llaftgcd w 


dth ar 




sources detected 


, as a hard X-ray point source ars 


i flagged wi 


th an 'x' 


, sources with an 


Optical AGN spectrum 


are flagged 


with an 


's', and sources with a UV quasar spcctru 


m are fit 


igged 


with 


a 'g'. 


































Optically unresolved sources /2 


< 


'.'6) with 


g — r 


< 0.1 


6 and u - 


- g < 1 


are classifi 






(QSO). 


Optically unresolved sources with u - 


- g > 1.75 or 


g — r 


< 0.6 and u - g 


> 


1 are 


classified 


as stars (STAR). Optically 


resolved sources (t- 


1/2 ^ o''6) are classified as gal; 


axles (GAL), 




cally resolved sources with a hard X- 


ray detection 


or an acti 


ve galaxy s; 


pectrum 



are classified as active galactic nuclei (AGN). 
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TABLE 9 

Optical Matches to UV Variable Sources in CFHTLS D4 



No. 


GALEX ID 


u 


9 


r 


i 


Z 




Sep 


flag'' 


class'' z 














(arcsec) 


(circsec) 






D4-1 


J221336. 1-180531 


19.459±0.000 


19.169±0.000 


19.208±0.000 


18.977±0.000 


19.001±0.001 


0.484 


0.3 


o 


QSO 


D4-2 


J221409.7-175652 


20.046±0.001 


19.894±0.001 


19.639±0.001 


19.541±0.001 


19.538±0.002 


0.465 


0.2 




QSO 


D4-3 


J221413.0-172909 


20.730±0.002 


20.584±0.002 


20.251d=0.002 


20.163±0.002 


20.072±0.004 


0.487 


1.2 


o 


QSO 


D4-4 


J221429.5-173007 


18.455±0.000 


18.126±0.000 


18.104±0.000 


18.126±0.000 


17.960±0.001 


0.459 


0.2 




QSO 


D4-5 


J221432.7-180158 


20.253±0.001 


19.812±0.001 


19.684±0.001 


19.638±0.001 


19.382±0.002 


0.484 


0.1 




QSO 


D4-6 


J221433.8-173007 


21.686±0.006 


20.844±0.002 


19.962±0.001 


19.524±0.001 


19.207±0.003 


0.831 


0.0 


o 


GAL 


D4-7 


J221508.9-175602 


25.067±0.075 


24.499±0.039 


24.584±0.048 


23.842±0.030 


23.438±0.069 


0.746 


0.8 


o 


GAL 


D4-8 


J221639.6-174427 


20.916±0.002 


20.396±0.001 


20.385±0.001 


19.995±0.001 


19.743±0.003 


0.532 


0.3 




QSO 


D4-9 


J221647.8-175454 


20.765±0.002 


20.358±0.001 


20.138±0.001 


19.673±0.001 


19.257±0.002 


0.595 


0.4 




QSO 


Sourct 


iS that are optically vari; 


able are flagged with s 


m 'o', sources detected as a hard X-ray 


point source are fl 




rees with 


an optical AGN specti 


:um are flagged 


with an 


'b', and sources with a UV quasar spectrum ar^ 


e flagged with a 'g'. 
















b Optic. 


ally unresolved sources ( 


^1/2 < ""^) ^^^^^ 9 - 


- r < 0.6 and u - g 


< 1 are classified 


as quasars (QSO). 


Optically unresolved 


sources w 


■ith -u — g > 


1.75 or g 


- r- < 0.6 and 


u - g > 


1 are classified as stars 


(STAR). Optically res 


olved sources (i-j^j 


> 0^'6) are classified as galaxies (GAL), and optically resolved sources with a hard X-ray 


detection or an 


active g£ 


ilaxy spectrum are elassi 


fied as active galactic 


nuclei (AGN). 

















TABLE 10 

Tidal Disri-ptiox Flare CAxnTnATES 



No. GALEX ID flare flag temp*' photo-2 spec-2 spec type*" 



Dl-1 


J022409.7-044007 


F 


o 


em 


0.52 - 


0.60 


0.436 


Sy2 


Dl-9 


J022517.0-043258 


T 


o 


Sb 


0.28 - 


0.36 


0.326 


early 


Dl-10 


J022518. 1-043155 


F 


o 


em 


0.48 - 


0.60 






Dl-14 


J022544.4-040102 


F 




em 


1.36 - 


1.48 






D2-9 


J100002.0+024216 


F 




em 


0.80 - 


0.88 


0.579 


Sy2 


D2-11 


J100014.9+022717 


F 




em 


0.00 - 


0.08 


0.727 


Syl 


D3-3 


J141721. 2+530916 


T 


o 


Sa 


0.48 - 


0.60 


0.355 


Sy2 


D3-8 


J141833.5+530733 


T 


o 


cm 


0.72 - 


0.80 


0.681 


? 


D3-13 


J141929.8+525206 


T 


o 


E/SO 


0.28 - 


0.36 


0.3698 


early 


D3-21 


J142107.7-I-530319 


F 


o 


em 


0.96 - 


1.04 


0.735 


Syl 


D4-7 


J221508.9-175602 


F 


o 


em 


0.96 - 


1.12 


0.189 


SN Iln 



^ Galaxy template from photo-z fitting routine. Emission line galaxies (either star forming or Seyfert 

galaxies) are labeled as 'em'. 

Spectra with no emission lines are classified as early tyjje galaxies (early). Spectra with diaf!,iiostic 
narrow-line ratios of ([O III] A5007)/(H/3) > 3 are classified as Seyfert 2 galaxies (Sy2) and spectra that 
show broad H/3 emission are classified as Seyfert 1 galaxies (Syl). 



TABLE 11 

Log of Optical Spectrsocopic Observations 



No. 


GALEX ID 


UT Date 


Telescope 


Instr. 


slit 


AA 


tcx 


z 


[O III]/H/3 


log(L([0 III])) 

oV \L J// 












(") 


(A) 


(sT 






(ergs s~^) 


Dl-1 


J022409. 1-044007 


2007-09-16 


Keck 


LRIS 


1.0 


5.7 


1800 


0.436 


4.0 


41.3 


Dl-9 


J022517.0-043258 


2006-09-26 


VLT 


FORSl 


0.7 


5.0 


2250 


0.326 




<37.8 


D2-9 


J100002.0+024216 


2007-05-18 


Keck 


LRIS 


0.7 


4.8 


2700 


0.579 


3.6 


41.1 


D2-11 


J100014.9+022717 


2007-05-18 


Keck 


LRIS 


0.7 


4.8 


1800 


0.727 


2.0 


41.4 


D3-3 


J141721. 2+530916 


2005-12-28 


MDM 2.4m 


CCDS 


2.0 


15.0 


2700 


0.355 


9.0 


40.4 


D3-8 


J141833.5+530733 


2007-04-19 


Keck 


LRIS 


1.0 


5.7 


2700 


0.681 


2.7 


40.4 


D3-13 


J141929.8+525206 


2007-04-19 


Keck 


LRIS 


1.0 


5.7 


2400 


0.3698 




< 39.6 


D3-21 


J132107.7+530319 


2007-04-19 


Keck 


LRIS 


1.0 


5.7 


3600 


0.735 


> 2.0 


40.6 



